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I.  INTRODUCTION 


During  the  grant  period  5/15/83  -  8/14/84,  the  efforts  of  our  group  were 
directed  towards  investigating  surface  scattering  of  ions  for  the  purpose  of 
obtaining  stimulated  emission  of  soft  x-rays,  as  well  as  scaling  of  free 
electroi  laser  to  x-ray  wavelengths,  and  nuclear  spectroscopy  in  the  x-ray 
regime. 

The  following  report  consists  of  three  sections  and  six  appendices.  The 
first  section  summarizes  scattering  of  ions  from  different  types  of  surfaces  and 
stimulated  emission  of  light  from  scattered  ions.  The  second  section  summarizes 
our  continuing  work  in  the  area  of  an  x-ray  free  electroi  laser.  Another 
continuing  effort  of  our  group  has  been  in  nuclear  spectroscopy  with  x-ray 
lasers  and  intense  fields.  The  work  in  this  area,  corresponding  to  the  relevant 
grant  period,  is  summarized  in  the  third  section.  The  papers  which  detail  the 
results  of  these  investigations  either  in  published  or  preprint  form  are 
reproduced  in  the  appendices. 


II.  SURFACE  SCATTERING  01’  IONS  AND  SOFT  X-RAY  EMISSION 


Coherent  excitation  of  ions  by  the  periodic  potential  of  a  crystal  proposed 

by  Okorokov  sometime  ago  arises  from  the  fact  that  an  ion  moving  in  a  crystal 

with  speed  v  experiences  a  time  dependent  potential  with  harmonics 

=  (2nvn)/d,  where  n  is  a  positive  integer  and  d  is  the  lattice  spacing  in  the 

direction  of  ionic  motion.  If  v  and  d  are  such  that  ft  matches  an  ionic 

n 

transition  frequency,  then  ions  are  resonantly  excited  by  the  n^-  harmonic.  The 
probability  that  ions  are  in  an  excited  state  becomes  a  periodic  function  of  the 
thickness  of  the  foil,  which  can  be  viewed  as  a  consequence  of  the  Rabi 
oscillations  of  ions  in  their  rest  frame.  One  can  expect  the  same  physical 
picture  to  hold  for  ions  scattering  from  crystal  surfaces  at  grazing  incidences. 
This  may  be  called  the  surface-Okorokov  effect  and  offers  interesting 
possibilities  for  obtaining  stimulated  emission  of  soft  x-rays  from  ions.  One 
can  coherently  excite  a  beam  of  ions  into  certain  high  lying  states  by  an 
adjustment  of  the  beam  velocity  to  the  crystal  periodicity,  and  obtain  an 
effective  population  inversion  relative  to  some  low  lying  states.  For  such  a 
purpose,  the  grazing  angle  surface  scattering  geometry  has  considerable 
advantages  over  the  beam  foil  geometry  in  that  the  damage  to  the  surface 
resulting  from  ion  bombardment  is  minimized  at  grazing  incidences  and  the  same 
surface  can  be  used  repeatedly.  Furthermore,  stopping  powers  of  surfaces  can  be 
lcwer  than  stopping  powers  of  bulk  materials,  which  means  that  ions  can 
resonantly  interact  with  surfaces  for  longer  distances  than  with  foils. 

We  analyzed  the  scattering  of  hydrogenic  ions,  such  as  Li+“,  Be*'*,  etc., 
from  crystal  surfaces  at  grazing  incidences,  and  showed  that  ions  execute 
extremely  rapid  Rabi  oscillations  in  their  rest  frame  (^1017  sec  1).  We  assumed 
that  the  ionic  beam  simultaneously  interacts  with  a  coherent  electromagnetic 
field  whose  frequency  matches  one  of  the  ionic  transition  frequencies 


corresponding  to  an  effective  population  inversion,  and  shewed  that  the 
inhomogeneous ly  broadened  gain  of  the  signal  can  be  substantial,  on  the  order  of 
1  cm"1,  for  reasonable  beam  currents  such  as  4.8  ampere/mm2 .  The  analysis  is  detailed 
in  the  preprint  in  Appendix  A,  which  is  accepted  for  publication  in  Phys.Rev.  Letters. 

In  the  free  electron  laser,  lasing  occurs  at  the  wavelength  X  =  2Xq/y, 
where  Xq  is  the  wiggler  wavelength  of  the  magnetic  selonoid  and  y  is  the  ratio 
of  the  electron's  energy  to  its  rest  mass  energy.  In  principle  one  can  decrease 
X  either  by  increasing  the  electronic  energy  or  by  decreasing  the  wiggler 
wavelength  Xq.  The  latter  method  is  not  very  practical  with  mechanical  systems 
used  in  free  electron  lasers.  There  are,  hewever,  systems  which  are  periodic  on 
atomic  scale,  and  which  might  be  used  as  short  wavelength  wigglers  for  electrons 
or  ions.  An  interesting  idea  is  to  use  anti ferromagnetic  crystals  as  wigglers 
for  ions.  These  crystals  have  strong  internal  magnetic  fields  that  are  periodic 
with  some  lattice  constant  and  can  be  used  to  excite  ions  coherently.  Excited 
ions  can  then  decay  radiatively  and  thus  convert  the  virtual  photons  of  the 
material  magnetic  field  to  real  photons. 

We  analyzed  the  problem  for  a  beam  of  hydrogenic  Be+3  ions  scattering  from 
the  surface  of  an  anti ferromagnetic  material  such  as  Fel’2  or  MnO.  In  the 
geometry  we  considered,  Be+3  ions  penetrate  just  below  the  surface  of  the 
material,  travel  some  distance  L,  and  then  come  out.  They  simultaneously 
interact  with  soft  X-Ray  signal  of  frequency  hws  ^  30  eV  which  matches  the 
frequency  ^32  of  the  n=3  -*■  n=2  transitions  of  Be+3.  We  first  considered  the 
magnetization  fields  inside  an  antiferromagnet  and  deduced  the  general  form  of 
the  fields  from  a  simple  argument.  We  then  obtained  a  simple  set  of  equations 
of  motion  for  the  ionic  amplitude  from  the  Schrodinger  equation  in  the  rest 
frame  of  an  ion.  These  equations  shew  that  ions  in  their  rest  frame  execute 
Rabi  oscillations  between  the  states  n  =  1  and  n  =  3,  if  the  ion  velocity  is 


properly  adjusted.  The  Rabi  frequency  fiR  is  proportional  to  Bq,  the  amplitude 
of  the  magnetization  of  one  of  the  sublattices  of  the  antiferromagnet,  and  to  v, 


the  ion  velocity.  f2R  is  on  the  order  of  4x10 1,1  sec  1  for  v  =  10  9  cm/sec,  which 
is  nearly  four  orders  of  magnitude  larger  than  the  spontaneous  radiative  decay 
rate  of  the  n  =  3  states  of  Be+3.  Because  ions  are  actually  moving  with 
velocity  v,  the  Rabi  oscillation  in  the  rest  frame  become  spatial  oscillations 
in  the  lab  frame.  If  LR  designates  the  spatial  interval  for  half  a  Rabi  cycle 
(i.e.,  Lr  is  the  minimum  distance  required  for  an  ion  to  travel  in  the  material 
in  order  to  get  the  ion  fully  excited  into  the  n  =  3  states,  if  it  is  initially 
in  the  ground  state),  then  LR  is  on  the  order  of  10~4  cm.  We  also  analyzed  the 
gain  of  the  signal  field  in  the  weak  signal  regime.  Due  to  the  coherent 
excitation  of  ions,  the  expression  for  the  gain  differs  from  the  usual 
small-signal  swept-gain,  and  have  a  number  of  novel  features.  It  turns  out  that 
only  for  a  narrow  region  of  detuning  between  u  and  . ^  can  one  get  finite  gain. 
The  maximum  gain  is  about  10  cm  1  if  the  ionic  density  is  on  the  order  of  10 1‘* 
cm  3.  Appendix  B  reproduces  a  preprint  which  gives  the  details  of  these 
investigations. 


Our  work  with  magnetic  surfaces  has  yielded  another  idea  which  might  lead 
to  a  new  class  of  lasers.  In  diatomic  molecules  formed  by  the  group  VI  elements 
(e.g.,  0.-,  Sp,  SO)  and  by  a  combination  of  group  V  and  group  VII  elements  (e.g., 
NF,  NCI,  PCI),  the  spin  selection  rule  AS  =  0  leads  to  highly  metatable  states 
which  can  be  very  useful  for  high  power  laser  systems.  A  well-known  case  is  the 


state  of  the  oxygen  molecule.  Due  to  spin  and  parity  symmetry,  its  decay  to 
the  ground  state  is  highly  forbidden  with  spontaneous  decay  rate  2.6  x  1(J~  ' 
sec  l.  This  fact  is  used  in  oxygen-iodine  lasers  to  store  the  energy  in  the  1  , 


state  of  0;,  which  is  then  collisionally  transferred  to  atomic  iodine  via  the 


reaction  0,  ( *  A^,)  +  I  (  P3/  ■> )  v  0,(  'y  )  +  I(rPj^,)*  Radiative  emission  takes 


place  on  the  iodine  through  its  decay  back  to  2P^?/  which  is  a  magnetic-dipole 


allowed  transition  with  a  spontaneous  transition  rate  of  7.7  sec 


-l 


The 


frequencies  of  3A  -  3E  and  2P,  ..  -  2P,,_  transitions  match,  which  makes  it 
possible  to  resonantly  transfer  energy  between  O2  and  I.  In  practice  it  is 
difficult  to  find  such  matching  partners  for  other  metastable  systems  which 
might  have  applications  to  high  pcwer  lasers.  For  this  reason,  as  well  as  to 
avoid  loses  involved  in  collisional  energy  transfers,  there  has  long  been  an 
interet  in  obtaining  radiative  emissions  directly  from  the  excited  metastable 
system  itself.  To  achieve  direct  radiative  emission,  it  is  necessary  to  break 
the  various  symmetries  of  the  metastable  system.  For  instance,  the  parity 
symmetry  can  be  broken  by  applying  static  electric  field  or  nonresonant 
radiation,  or  by  means  of  nonresonant  collisons  with  another  molecular  or  atomic 
system.  On  the  other  hand,  breaking  of  the  spin  symmetry  requires  application 
of  either  a  spatially  varying  magnetic  field  as  in  a  Stern-Gerlach  device,  or  a 
direct  spin-spin  coupling  another  system.  This  latter  possibility  can  be 
realized  by  putting  the  singlet  delta  oxygen  in  contact  with  a  magnetic  surface, 
where  the  electrons  of  the  singlet  delta  oxygen  will  couple  to  the  external 
spins  through  their  internal  magnetic  moments.  By  means  of  a  simple 
perturbation  analysis,  we  shewed  that  both  parity  and  spin  symmetries  are  broken 

and  that  an  electric  dipole  moment  is  induced  between  3A  and  3l"  which  is  on 

g  g 

the  order  of  10  2  Debye.  This  implies  a  spontaneous  transition  rate  of  about 
10 2  sec"1,  which  is  nearly  six  orders  of  magnitude  larger  than  the  corresponding 
rate  for  isolated  02.  The  analysis  should  also  hold  for  OF  which  has  an 
ultraviolet  transition  frequency  between  the  electronic  ground  states  and  the 


first  excited  singlet  delta  states.  This  work  is  already  published  (Phys.  Rev. 
A  29,  February,  1984).  The  paper  is  reproduced  in  Appendix  C. 


III.  FREE  ELECTRON  LASER  AT  X-RAY  WAVELENGTHS 

An  x-ray  free  electron  laser  must  operate  with  large  gain  per  pass  because 
of  the  poor  quality  of  available  mirrors  at  x-ray  frequencies.  Moreover,  the 
effects  of  inhomogeneous  broadening,  start-up  from  noise,  coherence  development, 
and  quantum  recoil  can  be  important  in  determining  whether  and  how  such  devices 
will  operate.  We  performed  one-dimensional  classical  and  semiclassical  analyses 
of  the  gain  regimes  for  x-ray  free  electron  lasers.  Based  on  these  analyses,  we 
determined  some  of  the  constraints  on  such  devices  imposed  by  effects  such  as 
diffractive  spreading  of  the  laser  beam,  energy  spread  and  emittance  of  the 
electron  beam,  and  transverse  variations  in  the  wiggler  field.  These 
constraints  are  discussed  in  the  preprint  reproduced  in  Appendix  D. 
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IV.  NUCLEAR  SPECTROSCOPY 

In  two  recent  papers  LPhys.  Rev.  C  TT_  (1984),  1229],  Reiss  claimed  that 
forbidden  nuclear  beta  decay  can  occur  in  the  presence  of  an  intense  but  low 
frequency  electromagnetic  field.  His  approach  treats  the  weak  coupling  to  first 
order,  and  the  electromagnetic  field  as  nearly  as  possible  to  all  orders.  Thus, 
he  describes  electron  states  by  Volkov  wavef unctions,  whereas  he  approximates 
nuclear  states  by  the  so-called  manentum-translation-approximation 
wavef unctions.  We  have  now  shown  by  explicit  calculation  that  within  this 
model,  the  total  beta  decay  is  essentially  independent  of  the  external  1 

electromagnetic  field.  This  conclusion  holds  both  xn  nonrelativistic  and 
relativistic  treatments  of  the  problem.  Moreover,  by  comparing  the  "r  .  2"  and 
"p  •  A"  form  of  the  interaction  Hamiltonian,  we  have  shown  that  the  wavefr  .on 
in  the  momentum  translation  approximation  is  the  unperturbed  state  in  the 
Coulomb  gauge.  This  work  is  reported  in  a  preprint  and  a  paper  which  appeared 
in  Phys.  Rev.  C  29,  March  1984.  These  are  reproduced  in  Appendices  E  and  F, 
respectively. 
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Abstract 


We  discuss  coherent  excitation  of  ions  incident  on  crystal  surfaces  at 
grazing  angles  and  show  that  substantial  gains  for  soft  x-rays  can  be  achieved 
for  reasonable  beam  currents. 


Coherent  excitation  of  ions  by  the  periodic  potential  of  a  crystal  was 

proposed  by  Okorokov  sometime  ago  [1]  and  has  been  observed  by  different  groups 

who  passed  highly  ionized  atoms  through  thin  foils  L2,3].  The  effect  arises 

from  the  fact  that  an  ion  moving  in  a  crystal  with  speed  v  experiences  a  time 

dependent  potential  with  harmonics  =  (2Trvn)/d,  where  n  is  a  positive  integer 

and  d  is  the  lattice  spacing  in  the  direction  of  ionic  motion.  If  v  and  d  are 

such  that  matches  an  ionic  transition  frequency,  then  ions  are  resonantly 
n 

excited  by  the  harmonic.  The  probability  that  ions  are  in  an  excited  state 
becomes  a  periodic  function  of  the  thickness  of  the  foil  [4],  which  can  be 
viewed  as  a  consequence  of  the  Rabi  oscillations  [5]  of  ions  in  their  rest 
frame.  When  one  goes  to  the  laboratory  frame  in  which  the  crystal  is  at  rest, 
the  Rabi  oscillations  in  time  become  thickness  dependent  oscillations  in  space. 
The  Okorokov  effect  depends  critically  on  the  finiteness  of  the  probability  that 
ions  keep  their  charge  states  constant  without  further  ionization  or 
neutralization  in  the  crystal. 

One  can  expect  the  same  physical  picture  to  hold  for  ions  scattering  from 
crystal  surfaces  at  arazing  incidences.  This  may  be  called  the  surface-Okorokov 
effect  and  offers  interesting  possibilities  for  obtaining  stimulated  emission  of 
soft  x-rays  from  ions.  One  can  coherently  excite  a  beam  of  ions  into  certain 
high  lying  states  by  an  adjustment  of  the  beam  velocity  to  the  crystal 
periodicity,  and  obtain  an  effective  population  inversion  relative  to  some  lew 
lying  states.  For  such  a  purpose,  the  grazing  angle  surface  scattering  geometry 
has  considerable  advantages  over  the  beam  foil  geometry.  The  damage  to  the 
surface  resulting  from  ion  bombardment  is  minimized  at  grazing  incidences  so 
that  the  same  surface  can  be  used  repeatedly  [6].  Furthermore,  stopping  powers 
of  surfaces  can  be  lower  than  stopping  powers  of  bulk  materials,  which  means 
that  ions  can  resonantly  interact  with  surfaces  for  longer  distances  than  with 


foils. 


In  this  paper  we  consider  the  scattering  of  hydrogenic  ions,  such  as  Li  +  2, 
Be+3,  etc.,  from  crystal  surfaces  at  grazing  incidences,  and  show  that  ions 
execute  extremely  rapid  Rabi  oscillations  in  their  rest  frame.  We  assume  that 
the  ionic  beam  simultaneously  interacts  ith  a  coherent  electromagnetic  field 
whose  frequency  matches  one  of  the  ionic  transition  frequencies  corresponding  to 
an  effective  population  inversion,  and  show  that  the  inhomogeneously  broadened 
gain  of  the  signal  can  be  substantial. 

The  geometry  of  the  problem  is  shown  in  Fig.  1.  An  ion  is  incident  from 
the  left,  moving  in  the  z'-direction  (primed  coordinates  refer  to  the  laboratory 
frame,  unprimed  ones  to  the  rest  frame  of  the  ion).  In  the  laboratory  frame, 
the  periodic  crystal  potential  can  be  described  by 

V^(x')  =  9(x0  -  x')Lq;  elG‘X'  (1) 

G  G 

Here,  Xo  designates  the  distance  between  the  c.o.m.  of  the  ion  and  the  surface 
(xo  is  negative  if  the  ion  is  above  the  surface;  x0  is  positive  if  the  ion 
penetrates  the  surface,  and  is  actually  inside  the  crystal).  G's  are  the 
reciprocal  lattice  vectors.  We  have  assumed  that  the  crystal  terminates  sharply 
on  a  plane  at  xq .  This  plane  need  not  coincide  with  the  plane  of  the  first 
layer  of  surface  atoms;  rather,  it  should  be  considered  as  an  effective  plane 
where  the  influence  of  surface  atoms  and  electrons  of  the  crystal  on  the  ion 
begins  to  be  appreciable. 

If  one  goes  to  the  ionic  rest  frame,  \T(x‘)  is  transformed  into  a  vector 
(A.)  and  a  scalar  potential  (V  ).  A,  is  purely  longitudinal,  and  can  be 
eliminated  by  a  gauge  transformation  on  the  wavef unction  of  the  ionic  electron, 
<f<(x,t)  =  exptiA(x,t)]  •  <j>(x,t),  such  that  VA  =  (e/hc)Ac<  One  finds  that 
satisfies  the  Schrodinger  equation 


ifi3t<)>  =  [(2m)  1  (p  -  eAsc  *)  -  Ze2jx|  1  +  +  V|j  ] <J> 


where 


Vi(x1)  =  ey9(x0  -  x)JsJ«l^  e1  X,  si  =  {G |  z  •  G  =  0}  , 

G 


V(1  (x,t)  =  eY_19(x0  -  x)Es||<t7  e1  X,  s||  =  {G |  z  •  2  /  0}  , 

G 


and  As  =  yAg  cos(ksz  -  wst)  is  the  signal  field,  y  is  the  relativistic  factor 

_j. 

(1  -  v2/c2)  2  and  essentially  equals  1  for  all  v's  of  interest  here.  is  a 

static  potential  that  depends  only  on  the  transverse  coordinates  of  the  ionic 

electron.  Note  that  going  from  the  laboratory  frame  to  the  ionic  rest  frame  is 

equivalent  to  taking  a  time  average  of  the  crystal  potential  along  the 

trajectory  of  the  ion;  Vj_  is  this  time  averaged  potential.  is  obtained  from 

this  perspective  in  Refs.  7  and  8,  and  the  level  shifts  induced  by  on  ionic 

spectra  in  foils  are  estimated  in  Ref.  8.  Here,  however,  we  will  ignore  V^, 

assuming  that  it  produces  negligible  perturbations.  For  our  discussion,  the 

important  term  is  the  time  dependent  potential  V||  .  If  the  basis  vectors  of  the 

reciprocal  lattice  space  are  G3,  22  and  S3,  that  is  G  =  IG\  +  m22  +  nG3,  where 

£,m,n  are  integers,  then  the  harmonics  of  Vy  are  given  by 
A  — ►  — ►  ■+ 

yvz  •  ( W»i  +  m3,?  +  nG3) .  In  the  following,  we  will  consider  the  case  in  which 
one  of  the  fundamental  frequencies  obtained  from  this  expression  equals  or  is 
near  an  ionic  transition  frequency.  Specifically,  we  will  assume  that  only 
n  =  1  and  n  =  3  levels  of  the  ion  are  involved  in  the  transitions  induced  by  the 
crystal,  and  that  an  effective  population  inversion  occurs  between  n  =  3  and 


n  =  2  levels  as  shown  in  Fig.  1.  One  can  then  ignore  all  G  s  m  the  set  sll 
except  for  those  which  give  a  constant  frequency  R  which  is  either  resonant  or 
nearly  so  with  the  ionic  transition  frequency  0131.  Thus  we  define  the  set  of 


vectors  g  such  that  g  =  xG^  +  yG^r  +  zYG^  and  |g*  £  I  =  constant  =  ft/v.  Ignoring 
the  signal  field  for  the  moment,  the  amplitude  equations  in  the  rotating  wave 
approximation  are 


i(fi-i03i)t 


IC100  =  e 


M 


100;32.m  35,m 


(3; 


■  i  (R-10  3 1 ) t 


ic,„  =  e 
3£m 


Ml00;3£m  C100* 


(31 


where 


MnS.m;n'  Z'm'  =  e(h'^ 


d*  un*m(x)0(x°  ‘  x)el8’X  Va'm'(x)  * 


and  u  ' s  are  the  hydrogenic  wave functions.  Assuming  c100(O)  =  1  and 

nS.m 


c  (0)  =  0,  the  solution  of  (3a,fc>)  is 
3«.m 


ClOO 


=  (flV1*1  *  -  R~  e"in  t]/(n+  -  R")  , 


(ein+t  -  eiart)/[Q*  - 


c _ .  (t)  =  M.  .  Ie‘**  t  -  e1^  *  /(R  -  R  )  , 


(3c 


(4 


(41 


where  jr(  ta)  =  U31  "  n)/2  ±  C(u31-  Q)1/*  +  2  ^  nRJlm  “  lMi00;3(Un'  is  thi 

Rabi  frequency. 


If  the  ion  penetrates  the  surface  and  |x0l  ^  1A  >  (a^/Z) ,  one  can 
essentially  ignore  the  step  function  in  (3c)  since  u  *s  316  localized  to  a 
radius  of  approximately  (na  )/Z.  One  then  finds 

D 


100 : 300 


=  3/37  e(4hY)_1L qi  8W  +  3)  (B*  +  1)‘ 


10  0 ; 31m 


=  i/7  e(2hY/3)_1^33(e2  +  2)  CS|  +  D"3Ylm(g)  , 


100;51m  2hY 


Lit.  ¥2 


-  (*)6g 


6  (363  +  4)  pr  rfs  +  f]  rfs  -f 

*  2-& -4-  -  r  — r^Yr-r 

6 (62  +  l)4  y  s-0  27  s+1  r  s+4 


•*! 


■n  *  ‘j 


where  6  =  4Z/(3ga  )  and  the  superscript  p  is  meant  to  imply  a  penetrating  ion. 

g  B 

Note  that  only  m[oo-3io  describes  30  electric  dipole  transition;  the  other 
transitions  violate  the  dipole  selection  rules  since  they  arise  from  higher 
order  multiple  transitions.  This  is  in  contrast  to  the  formulation  of  the 
Okorokov  effect  in  Rfs.  4  and  9,  where  the  dipole  approximation  is  assumed  to 
hold.  In  justification  of  (5a-c),  we  note  that  whether  one  can  make  the  dipole 
approximation  depends  on  the  ratio  of  the  crystal  period  to  the  radius  of  the 
ion's  electronic  orbital,  that  is  on  whether  |g|(a„n/Z)  <<  1.  For  typical 
crystals  one  can  easily  obtain  |g|(a  n/Z)  -v  1  even  for  low  lying  orbitals.  The 
dipole  approximation  breaks  down.  Incidentally,  this  fact  can  be  exploited  to 
populate  those  states  which  are  not  dipole  connected  to  the  ground  state  and 
therefore,  relatively  long-lived.  After  leaving  the  surface,  the  ion  remains  in 


a  superposition  of  these  long-lived  excited  states  with  the  ground  state.  This 
superposition  is  stable  relative  to  radiative  decays  to  the  ground  state. 

If  one  takes  the  pseudopotential  values  of  <y  that  are  used  in  structure 

I 

calculations  LlOj,  and  uses  just  the  smallest  and  the  next  smallest  g's  in  the 

sums  of  (5),  one  finds  for  a  crystal  like  Si  that  ,nJ»  IM^  __  I  -v  1016 

sec"1  and  lM[oo-3lo^  ^  ^1?  sec 1  •  The  resonance  condition  vz  •  (Gi  +  G2  +  G3) 

^  “3i  ^  193  eV  can  be  satisfied  for  a  Be+3  ion  moving  with  speed  v  n,  109 

cm/ sec.  It  follows  that  the  ion  can  be  completely  excited  after  traveling  a 

"Rabi  distance"  £  =  ttv/Q  v  3  -  30A. 

R2m  R£m 

An  interesting  feature  of  (5b,c)  is  that  the  contributions  of  different 

components  are  proportional  to  Y  (g).  There  are,  therefore,  orientation 

£m 

dependent  interferences  among  cp  components.  By  varying  the  direction  of  the 

I 

ionic  beam,  it  is  possible  to  pick  different  sets  of  g-vectors,  and  thus  obtain 

information  about  different  <V^  components  from  a  comparison  of  Rabi  frequencies 

g 

corresponding  to  these  sets. 

We  now  consider  the  coupling  of  the  ion  to  the  coherent  signal  field  and 

treat  this  coupling  as  a  perturbation  imposed  upon  the  evolution  described  by 

(4a, b).  In  the  calculation  of  the  gain  the  steps  are:  a)  determine  c2im(t)  to 

first  order  in  A  using  c_„  (t)  given  by  (4b);  b)  calculate  the  transverse 

current  density  resulting  from  (3£.m)-*(2e.,n' )  transitions;  c)  substitute  this 

transverse  current  density  into  the  f-axwell  equation  for  Ag  to  obtain  a 

dispersion  relation  between  a  real  us  and  a  complex  k^ .  The  imaginary  part  of 

ks  gives  the  small  signal  gain  for  ions  moving  at  a  fixed  velocity  v.  Finally, 

d)  multiply  the  resulting  gain  with  the  probability  distribution  of  ion 

velocities,  p.  (v),  and  replace  n  by  v|z*g|.  We  assume  that  the  width  of 
ion 

pion(v)  is  sufficiently  narrow  as  not  to  effect  the  overall  choice  of  the  set  of 
g's,  and  integrate  over  the  interval  -  »<  v  <  ®.  in  the  limit  that  the 


lifetimes  of  the  states  in  the  n  =  2  shell  go  to  the  resulting 
inhomogeneous ly  broadened  gain  is  given  by 


e  j2  1>t  V°ion<V> 

“  “s”2c  *“°  U“s  ‘  "32)2  *  “k1 


where  N.  is  the  density -of  ions,  P  is  the  probability  that  the  charge  state 
ion  cs 

of  ions  remains  unchanged,  and  v^m  is  a  tuning  velocity  given  by 


-  u>  -  U)__  -  -  +  w_. 

»  ->i  S  32  0)  oj,-  31 

2  •  g|  s  J 


D  is  a  coefficient  related  to  the  momentum  matrix  elements. 


n  _  V  Y  <?  <?*  M  M* 

“  LV m'  2rm' ;3£"m"  2J,,m*  ;3Hm  100;3£m  100; 3«.' 


where 


^24m; 3£'m'  =  u2l»(-i7)u3£’m' 


=  -(Z/aB)5m,0(5mM  +  +  ,2^  ) 


Here  m  =  213  •  32  •5"10=  7.6x10  "3  and  n2  =  21B  •  33  •  5~10  =  0.72  are  pure 
numbers  [11 J.  A  representative  gain  curve  vs.  ws  is  depicted  in  Fig.  2  just  for 
one  particular  (icm)  in  the  sum  (6a),  and  for  a  gaussian  probability  distribution 
of  average  velocity  v0  and  width  Av:  Pion(v)  =  (Av/n)  1  expL-(v-vo)  2/Av2J.  As 
seen  from  the  figure,  the  gain  curve  has  double  peaks  surrounding  the  exact 
resonance  at  «  =  uj,~,  where  it  vanishes.  Q..  is  the  sum  of  several  such  curves 


A.  9 

Which  have  differing  heights  and  widths,  but  always  vanish  on  exact  resonance, 

“s  =  w32* 

Another  interesting  feature  of  £  given  by  (6a)  is  that  instead  of  the 

usual  population  inversion  factor  (N  -  N,  ),  it  has  an  effective 

upper  lower 

population  inversion  factor  Nionr2Rim/C(^s  -  ^3->)2  +  Oj^].  Thus  all  ions 

participate  in  the  photon  emission  and  absorption  processes.  The  degree  of  this 
participation  is  determined  by  the  coherence  of  ions  and  the  corresponding  Rabi 
frequency  caused  by  the  crystal. 

Going  back  to  the  Be+3  example,  one  finds  (  assuming  <v  .5  ) 

Ss  'v  (Nion/1012  cm  3 )*(v/av)x10~4  cm'1  for  a  signal  of  ^  ^  30  eV.  Thus,  for 
Av/v  n.  10~4  and  an  ionic  current  density  j.^  %  4.8  ampere/mm2,  one  has 

G$  v  1  cm'  1 . 

As  an  ion  moves  inside  the  crystal  or  on  its  surface,  it  loses  energy  to 

the  material  electrons,  and  its  c.o.m.  slows  down.  This  can  cause  the  ion  to 

move  away  from  the  resonance  with  the  crystal  field.  Scaling  from  the  stopping 

powers  of  protons  in  bulk  materials  [12],  one  finds  S  %  109  eV/cm  for  Be+5  with 

v  'v,  109  cm/sec  for  a  Si  target.  For  grazing  angle  scatterings  one  can  expect  S 

to  be  lower  since  the  number  of  material  electrons  may  be  depleted  near  the 

effective  plane  at  xq  due  to  the  formation  of  inversion  layers  further  down  in 

the  interior  [13 ].  It  is  also  possible  to  choose  the  beam  direction  in  such  a 

way  that  ions  can  move  between  the  crystal  planes,  avoiding  heavy  concentration 

of  material  electrons.  Such  a  channeling  of  ions  can  decrease  S  by  orders  of 

magnitude  [14].  For  an  ion  moving  the  Rabi  distance  z  ,  the  minimum  energy  loss 

is  E,  v  Se.D .  Taking  S  -v  103  eV/cm  for  channeled  ions  and  0^3-  30A,  one 

finds  AE.  /E.  v  10  6  -  10  5  for  Be*0,  which  is  quite  small, 
loss  ion 
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FIGURE 

FIG. 1. 
FIG. 2. 


CAPTIONS 


The  schematic  of  the  grazing  incidence  surface  scattering. 
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For  illustrative 


The  gain  for  a  particular  (jim)  vs.  5  =  (u  -w.,-)  (|z*g|iv) 

o  S  jZ 

purposes,  |  z*g[  Av) -1  and  oc=(Ui31~v0 1  z*g [ )  ( | z*g  |  av) -1  are  taken 


to  be  3  and  1, respectively. 
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abstract 

We  show  that  hydrogenic  ions  scattering  from  antiferromagnetic  surfaces  execute 
rapid  Rabi  oscillations  in  their  rest  frames,  which  can  be  used  to  amplify  X-rays 
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I.  INTRODUCTION 

The  subject  of  the  present  paper  was  inspired  by  the  free 
electron  lasers.  In  free  electron  lasers,  the  laser  action,  or 
amplification,  occurs  at  the  wavelength  X  =  2X^/y  ,  where  X^  is  the 
wiggler  wavelength  of  the  magnetic  selonoid  and  Y  is  the  ratio  of  the 
electron's  energy  to  its  rest  mass  energy  [1],  In  principle,  one  can 
decrease  X  either  by  increasing  the  electronic  energy,  or  by 
decreasing  the  wiggler  wavelength  i  .  The  latter  is  not  very 
practical  with  the  mechanical  systems  that  are  used  in  free  electron 
lasers.  There  are,  however,  systems  which  are  periodic  on  atomic 
scale,  and  which  might  be  U3ed  as  short  wavelength  wiggler s  for 
electrons  or  ions.  A  particularly  interesting  idea  is  to  U3e 
antiferromagnetic  crystals  as  wigglers  for  ions.  These  crystals  have 
strong  internal  magnetic  fields  that  are  periodic  with  some  lattice 
periodicity  [2-4]  and  can  be  used  to  excite  ions  coherently.  Excited 
ions  can  then  decay  radiatively  and  thus  convert  the  virtual  photons 
of  the  material  magnetic  field  to  real  photons,  just  as  in  the  case  of 
free  electron  lasers. 

In  the  following  we  discuss  the  physics  of  this  problem  in  detail 
for  hydrogenic  ions.  Specifically,  we  consider  a  beam  of  Be*3  ions 
scattering  from  the  surface  of  an  antiferromagnetic  material  like  FeF 
or  MnO.  Be*3  ions  penetrate  just  below  the  surface  of  the  material, 
travel  for  some  distance  L,  then  come  out.  They  simultaneously 
interact  with  a  soft  X-Ray  signal  of  frequency  hu^  ^ 30  eV  which 
matches  the  frequency  ^32  of  the  n-3  -*•  n-2  transitions  of  Be*3.  In 
Section  II,  we  first  consider  the  magnetization  fields  inside  an 
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antif erromagnet  and  deduce  the  general  form  of  the  fields  from  a 
simple  argument.  We  then  obtain  a  simple  set  of  equations  of  motion 
for  the  ionic  amplitudes  from  the  Schrodinger  equation  in  the  rest 
frame  of  an  ion.  These  equations  show  that  ions  in  their  rest  frame 
execute  Rabi  oscillations  [5]  between  the  states  n=l  and  n=3 ,  if  the 
ion  velocity  is  properly  adjusted.  The  Rabi  frequency  ft  is 

K 

proportional  to  BQ,  the  amplitude  of  the  magnetization  of  one  of  the 
sublattices  of  the  antif erromagnet ,  and  to  v,  the  ion  velocity.  is 

on  the  order  of  4x10 13  sec  1  for  v»109  cm/sec,  which  is  nearly  four 
orders  of  magnitude  larger  than  the  spontaneous  radiative  decay  rate 
of  the  n-3  states  of  Be+3[6].  Since  ions  are  actually  moving  with 
velocity  v,  the  Rabi  oscillations  in  the  rest  frame  become  spatial 

oscillations  in  the  lab  frame  [7].  If  L  designates  the  spatial 

R 

interval  for  half  a  Rabi  cycle  (i.e.,  L  is  minimum  distance  required 

R 

for  an  ion  to  travel  in  the  material  in  order  to  get  the  ion  fully 
excited  into  the  n=3  states,  if  it  is  initially  i3  in  the  ground 
state),  then  1^  is  on  the  order  of  10' 4  cm.  We  also  give  a  detailed 
discussion  of  the  validity  of  the  model.  In  Section  III  we  calculate 
the  gain  of  the  signal  field  in  the  weak  signal  regime.  Due  to  the 
coherent  excitation  of  ions,  the  expression  for  the  gain  differs  from 
the  usual  small- signal  swept-gain  [8],  and  has  a  number  of  novel 
features.  It  turns  out  that  only  for  a  very  narrow  region  of  detuning 
between  and  uj32  can  one  get  finite  gain.  The  maximum  gain  is  about 

10  cm  ’  1  if  the  ionic  density  is  on  the  order  of  1014  crn-"^.  Section  IV 


gives  a  further  discussion  of  the  results  and  our  conclusions. 
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II.  ANT I FERROMAGNETIC  SURFACE  SCATTERING 

The  schematic  of  ion  scattering  from  an  antiferromagnetic 
material  is  shown  on  Fig.  1.  A  beam  of  Be+3  ions  grazes  the  surface 
of  an  antiferromagnetic  material  at  an  angle  typically  less  than 
0.5°  [9],  penetrates  the  surface,  and,  after  traversing  the  material 
for  some  distance  L,  comes  out.  Simultaneously,  the  beam  interacts 
with  a  soft  X-ray  signal  of  frequency  hws  ^  30  eV,  which  corresponds 
to  the  transition  frequency  between  the  n=3  and  n=*2  levels  of  Be  +  3. 
This  is  the  signal  that  we  want  to  amplify.  The  signal  propagates  in 
the  same  direction  as  the  ion  beam.  In  thi3  section  we  analyze  this 
problem  with  the  aid  of  a  simple  model.  We  first  discuss  the 
description  of  the  fields  of  the  antiferromagnetic  material,  then 
write  a  simple  set  of  equations  of  motion  for  Be+3  ion  and  explore  the 
implications  of  these  equations. 

An  antiferromagnetic  crystal  is  composed  of  two  interpenetrating 
identical  sublattices  which  are  shifted  relative  to  each  other  by  some 
spacing  d.  In  each  sublattice,  spins  are  localized  at  the  lattice 
points  and  are  parallel  to  each  other;  however,  they  are  antiparallel 
to  the  spins  of  the  other  sublattice.  The  vector  potential  induced  by 
the  intrinsic  magnetic  moments  of  these  spins  has,  therefore,  the 
periodicity  of  the  sublattice  and  can  be  written  as 


(x')=Z 

-+* 

G 


■  x ' 


X(G) 


(2.1) 


» 


» 


where  G  refers  to  the  reciprocal  lattice  vectors  of  the  sublattice. 

The  primes  are  meant  to  indicate  that  these  quantities  are  in  the  rest 
frame  of  the  material  medium.  The  presence  of  a  surface  breaks  the 
translational  symmetry  normal  to  the  surface  and  modifies  (2.1). 
However,  this  modification  is  relatively  simple  once  one  has  an 
expression  for  A (G) .  To  see  what  the  form  of  A(G)  is,  we  use  the 
simplified  picture  of  rigidly  fixed  spins.  Each  localized  spin  then 
gives  rise  to  a  magnetic  dipole  field,  the  superposition  of  which 
yields  the  overall  vector  potential: 


(x,)=gsyB5:- 


U* (x1 -R) 
1  X •  -Ri  3 


U* (x 1 -R-d) . 
■  — v  -►  -*■  -  3  J 
|  x ' -R-d | J 


(2.2) 


where  gsis  the  gyromagnetic  ratio  of  a  given  spin,  uB  is  the  Bohr 
magneton,  and 

R=£a+mb+nc  ,  l , m , n=integers  , 


are  the  lattice  vectors  of  one  of  the  sublattices.  u  designates  the 
direction  of  magnetization  of  this  sublattice.  The  Fourier  transform 


(2.4) 


= ( V/N ) ( e  1^‘^)  X(q)  , 

R 

where  V  is  the  volume,  N  is  the  number  of  the  sublattice  points,  and 
A(q)  is  defined  by 


£(q)=2TTN(l-e-1^) 


(v|q|2)  1{e+[(c_*u)  ( z *q) - ( z • y )  (e_*q) ] 


A  A  A  A  A  A  A  ^ 

-e_[ (e+*u) ( z  •  q)  -  ( z  •  y ) (e+*q) ] 

A  A  A  A  A  A  A 

+  z[ (e  +  *y)  (e_*q) - (e_*y)  (e  +  »q) J }  , 

(2.5) 


and 

e±=x' ±iy '  .  (2-6) 

When  the  number  of  sublattice  points  N  is  sufficiently  large  (N-«o), 
the  sum  over  R  gives 


Q-+ 


Z  e“iq’R  =NZ  5+ 


Inverting  the  Fourier  transform,  we  then  find 


(2.7) 


q 

-  Z,  ei2'J'  *<5)  . 


(2.8) 


Consider  now  the  special  case  in  which  u  is  along  x,  the  reciprocal 

/v 

lattice  is  cubic,  2  coincides  with  one  of  the  lattice  directions,  and 
is  midway  along  R3  such  that 


„  ,2tt,  ...  .  Trn 

G3=  ( 2^)  x  (integer) =  g— . 


(2-9) 


The  factor 


1  -  e"lG’d  =  1  -  (-i)n 


(2.10) 


is  nonzero  only  for  odd  n,  and  (2.8)  becomes 


Am*  (x')=  16::gsuBNV  e 

G 


K 


x  w  -*• 


n 


(K2+G2) 

n 


^  G,cosK  z' 

sinK  z'  -iz' - 2 —  ]  , 

n  2-2 

<VG  > 


(2.11) 


where 


K  =  ^.(2n+l) 
n  a 


(2.12) 


+ 


y ' 


m,£=integers . 


(2.13) 


Let  the  ion  beam  propagate  in  the  z-direction  with  a  uniform  cross 
section.  A  given  ion  then  samples  a  vector  potential  that  is 
approximately  the  average  of  (2.11)  over  the  transverse  coordinates. 
When  such  a  spatial  average  is  taken,  only  the  -  0  terms  survive  in 
(2.11),  and  the  vector  potential  becomes 


V<1‘,)~y’B0  Co 


sink  z ' 
n 


(2.14) 


where  B  3 


16-Ng  ..  /v.  The  corresponding  magnetic  field  is  given  by 
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I 


n=0 


cosK  z' 
n 


(2.15) 


When  a  surface  is  introduced,  the  basic  modification  occurs  in 

(2.7) .  Let  the  surface  be  the  y-z  plane.  We  separate  the  sum  in 

(2.7)  into  two  parts,  one  part  is  a  sum  over  the  sublattice  points  on 
the  y-z  plane  and  the  other  part  is  a  sum  over  the  sublattice  points 

A 

along  the  x-direction: 


r  e-i5*  ^  -ty  -i (mbq  +2ndg.)  -iq at 

V  "(Zm,ne  3)lZte  1  ) 


(NSZG  G  *q,,G ,fiq-,G  }  pl(c5^ 


(2.16) 


where  Ns  is  the  total  number  of  sublattice  points  on  the  y-z  plane  and 
pjtqi)  is  the  spin  line-density  fluctuation  normal  to  the  surface. 

Note  that  N  pj(0)  is  the  total  number  of  sublattice  points.  Eq.  (2.8) 
is  now  replaced  by 


i(qix’+G  y’+G  z')0  (a.) 

[  — - — 1 

l  ^  i  rt\  J 


■P1  (0) 


^ (q1,G2,G3) 


(2.17) 


If  one  ignores  the  spatial  variation  in  transverse  coordinates  and 
takes  an  average  over  them,  (2.17)  reduces  to  (2.14).  If  the  beam  is 

A 

of  small  cross-section  or  of  small  thickness  in  the  x-direction,  then 


transverse  variations  are  important  and  (2.8)  and  (2.17)  give 
different  results.  In  the  following  we  ignore  the  transverse 
variations  in  the  vector  potential  and  use  (2.14).  This  simplifies 
the  algebra  and  allows  us  to  concentrate  on  the  essential  physics  of 
coherent  excitation  of  ions. 

In  order  to  couple  the  electron  of  Be+3  to  the  antiferromagnetic 
vector  potential,  one  needs  to  transform  the  field  given  by  (2.14)  to 
the  rest  frame  of  an  ion  moving  with  velocity  v  in  the  z-direction. 

In  the  rest  frame  of  the  ion, 

A^x^-yyBq  i"=0  kn_1  sin(knz+u>nt)  ,  (2.18 

where 

Y=(1-v2/c2)'1/2  , 

*cn=YKn ' 

a)  =vvK 
n  n 

Thus  the  ion  experiences  a  time-dependent  field  in  its  rest  frame, 
whose  harmonic  frequencies  can  be  matched  to  one  of  the  ionic 
transition  frequencies  by  varying  the  beam  velocity.  In  particular, 
if  one  matches  one  of  w  to  <j31,  the  transition  frequency  between  the 
ground  state  and  the  n-3  states,  then  the  magnetic  field  of  the 
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material  can  resonantly  pump  the  ions  to  the  n-3  states.  The  ion  then 
executes  Rabi  oscillations  in  its  rest  frame.  The  frequency  matching 
wn  ■  “3i  requires  that 


“31  “31 


For  Be+ 3 ,  u31  ~  3xl017sec'1.  For  an  ant if erromagnet  like  MnO,  d  ~  4A. 
Substituting  these  values  into  (2.19),  one  finds  v(n-O)  %  3x10 9  cm/sec 
and  v(n»l)  *  109  cm/sec,  etc.  Finally,  if  the  matching  conditions  are 
met  for  the  beam  and  ions  are  coherently  excited  to  the  n-3  states, 
then  some  of  these  excited  ions  can  decay  into  the  n-2  states  by 
emitting  hu32  ~  30  eV  photons,  which  will  amplify  the  signal.  It  is 
clear  that  the  proposed  scheme  uses  the  same  type  of  excitation 
process  as  in  the  Okorokov  effect  [7,10-14]. 

In  order  to  simplify  the  equations  of  motion  we  assume  that 
pumping  is  strong  but  that  the  signal  is  weak.  We  first  consider  the 
Rabi  equations  for  the  n«l  and  n-3  states  coupled  by  the  material 
field  alone.  We  then  use  the  solutions  of  these  equations  to 
calculate  the  amplitudes  of  the  n-2  states  which  are  coupled  to  the 
n-3  states  by  the  weak  signal.  Let  the  ion  wave  function  in  its  rest 
frame  be 

-+  -iE  t 

•^  ( x ,  t )  =  y  C  (t)e  n  u  (  x  J  ,  (2.20) 

nim  n *.m 
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t 


where  u  ,  (xj  are  the  hydrogenic  wavef unctions  for  Be*3,  We  assume 
that  the  frequency  matching  condition  is  satisfied  only  for  the  nth 
harmonic  of  the  material  field.  Neglecting  the  lifetimes  and  the 
counter-rotating  terms,  one  then  has  the  Rabi  equations  for  n=l  and 

■v  3 

n=3  states  of  Be  : 


d 

dT  Ci°° 


i  e  y  B  r\ 
me  k 

n 


)  t 
n 


I 

£ '  m  * 


P  ( 1 00 ; 3  £  *  m ' 

y 


k 

n 


c3£  •  m* 


(2.21) 


) 


i  e  y  B  ^ 
me  k 

n 


i (aj l-M  }t 

e  n  Py(5zm; lOOj -knz)cl00  , 


(2.22) 


a 

a 


where 


P (n£m ; n ' £ ’ m '  k) 


d  x  e 


-  i  k  *  x 


n-£m 


('iV)un-£'m- 


(2.23) 


We  have  kept  just  the  p.A  terms  in  the  above  equations.  Performing 
the  indicated  integral  in  (2.23)  for  n=3  and  n=l  states,  one  obtains 
the  matrix  elements 


P  (3im; 100  |kr) 


(5 


m ,  1 


1F1 


(O 


+  i  6  -v  F  .,  ( £  )  ] 

(2.24) 
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where  Z  is  the  nuclear  charge  of  the  hydrogenic  ion  and  a  is  the  Bohr 

B 

radius.  The  other  quantities  are  defined  as  follows: 


1 6  Z  2  ♦  9  k2  a  2 

D 


(2.25) 


Fj  (5) 


^  (3+4;2)  +  c  3  (  t  _  r  2  1  l  m-j)  U-O* 

yrr  ^  5/-  *  =  °  r  ( i+4) r  C  i  + 1 )  2 1 


(l  +  O 


(2.26) 


F2CO  =  i5(l-0)’i[(6-^)0  *  (^-1)] 

4  0/2 


(2.2 7) 


Assuming  that  the  ions  are  initially  all  in  the  ground  state  (i.e.r 

c,,„  -  1,  c  -  0  at  t  -  0),  the  solutions  of  (2.21)  and  (2.22)  are 

1  J  3  1m 


c  1  0  0  (  1  )  " 


n-  -in  t , 

Q  e  )  , 


(2.28) 


(r.  -z  ) 


e>  B .  ~f‘^+  • o“ f 

c_,  ft)  - - : -  P  (31m;  1  0  0  |  -  k  z)  e1"  1  -  e1"  * 

jlm  ,  .  „  -  ,  -  »  v  n  I 

mck  (Z  -  Z  )  v 

n 


(2.29) 


where 


OJ  3  1  -  01 


w  3  1  " 


(2.30) 


Qr  is  the  Rabi  frequency  given  by 


e  Y  B  0  D  q  e  B  o  D  0  y  v 


(2.31) 


if  tCFl(€n))2  +  (F2Un))2,-*  , 

D 


(2.32) 


When  w31  =  un  exactly,  the  populations  of  upper  and  lower  states 
oscillate  with  the  Rabi  frequency  fiR.nR  depends  linearly  on  the 
amplitude  Bo  of  the  magnetic  induction  of  the  material.  For 
antif erromagnets  FeF;  and  MnO,  B0  ^’2x  107G  and  5.4xl06G,  respectively 
[15].  Taking  Bg  -  5xl0®G  as  typical,  and  assuming  that  the  zeroth 
harmonic  is  matched  to  u>  3  j  ,  we  find  D0  ^9x10 7  cm"1  and 

—  8.9xl013  sec  1  for  2=4.  Thus  if  the  ion  moves  for  a  distance  of 
LR  =  ttv/3  r  in  the  material  with  constant  velocity  v  -  3x10 9  cm/sec, 
then  the  population  of  ions,  which  is  initially  in  the  ground  state, 
i3  completely  excited  to  the  n=3  states.  Table  I  lists  some  of  the 
relevant  numbers  for  matching  ^31  to  some  of  the  harmonics. 

It  follows  from  (2.24)  that  only  certain  states  of  the  n=3  shell 


are  coupled  to  the  ground  state.  These  are  the  states  with  £=1  and 
£=2,  and  m=±l.  The  coupling  of  <-=2  states  obviously  arises  from 
higher  multipole  transitions  than  the  dipole,  since  the  change  in  the 


angular  momentum  A  £-2.  However,  F  2 ( - )  which  corresponds  to  Af  =  l  in 

(2.24)  also  has  contributions  from  the  higher  multipoles.  To  be 

precise,  the  electric  dipole  transitions  correspond  to  the  limit  ^  1 

of  P  .  In  this  limit  the  sum  in  (2-26),  as  well  as  ?2,  vanish,  and 

the  first  group  of  terms  in  (2-26)  yields  Fj(l)  =  0.1.  The  importance 

of  the  contribution  of  the  higher  multipoles  depends  sensitively  upon 

the  size  of  k  via  the  relation  (2-25). 
n 

Let  us  now  consider  the  amplitudes  of  the  n=2  states.  Let  the 
signal  field  be  described  by 


A.  =  yA  cos  (k_  z  -  u  t)  . 


(2.33) 


If  the  lifetime  of  the  n-2  states  is  r2»  then 


— -  +  r  7  c  =  - 

dt  z  2  £m  me 


leA  -i(co32-u)  )t 
s  J  ^  s 

-  e 


(2.34) 
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Substituting  the  solution  for  c3^  ,  from  (2.29)  and  integrating  the 

result  (assuming  c . >  -0  at  t=0),  one  finds 


e  <  B  n  A  _ 


C2£m(t)  = 


— - I -  [I£,m,Pv(^m;5^m'  |-k  z)P  (3Z'm'  ;100)-k„;)] 

m2  c  2  (Q-n)k  y  sy 
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P  i  ( P.  -'jJ32+^s)t  i  (Q  -OJ32+U  )  t 

e  e  s 


+  r  9  -  p 

|_  n  -(1)32+111^-1  — 2  -0)32+('J  -i 


2  -i 


(2.35) 


We  note  that  k  is  at  least  two  orders  of  magnitude  less  than  k  ,  and 

s  n 

therefore,  one  can  let  ksafi  -*■  0  and  use  the  dipole  approximation  in 

the  evaluation  of  the  momentum  matrix  element  with  k  .  In  thi3  limit, 
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P  ( 5  £m ; 2  £ ' m 1 
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|0) 


0(6m,  1 
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m ,  -  1 
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[5£,l6£',0/nl 


,26£M^2] 


(2.36) 


where  m  and  n2are  pure  numbers  given  by 


B.  16 


n 


1 


2  1  3  „2 

■—  1 0 °  =  7. 55x10" J  , 


.,18  ,3 

n  2  =  "ThT-  *  °-72* 


ft 

(2.37) 


It  is  worth  pointing  out  that  the  solution  (2.35)  indicates  a  Raman 
type  coupling,  since  the  amplitude  is  proportional  to  A^A  as  in  the 
Raman  scattering. 

In  the  preceding  discussion  we  assumed  that  the  velocity  of  the 
ion  remains  the  same  as  at  the  time  of  the  entry  of  the  ion  into  the 
material.  This  is  of  course  a  gross  approximation.  As  soon  as  the 
ion  enters  the  surface  region,  it  begins  to  lose  energy,  particularly 
to  the  electrons  of  the  target  material,  and  slows  down.  We  can 
estimate  the  error  involved  in  the  assumption  of  constant  ion  velocity 
from  the  stopping  power  of  electrons  in  solids,  which  depends  on  the 
speed  of  the  ion,  its  charge  number  Z,  and  the  charge  number  of  the 
target  atoms  ZT  [16-18].  The  stopping  power  scales  as 


S(Z, 


(2.38) 


where  is  the  stopping  power  when  the  projectile  is  a  proton,  and 
Z*(v)  and  Zp(v)  are  the  effective  charges  of  the  ion  and  the  proton, 
respectively,  in  the  material  [16] .  These  effective  charges  arise 
from  the  screening  effects  due  to  the  electrons  of  the  material  and 


ft 


> 


are  dependent  on  the  projectile  speed  [20-22].  Let  the  target  atoms 

2 

be  Mn  and  v  -  10 9  cm/sec.  For  this  speed  S  *  10"14eV.cm  /atom, 

P 

Z*/(ZZp)  **  0.6,  and  therefore,  S  ^  6x10  14  eV.cm2/atom.  The  number  of 
target  atoms  per  unit  volume  in  MnO  is  on  the  order  of  1022  cm’3, 

which  translates  into  an  energy  loss  per  unit  length  N^S  %  9x10 8 
eV/cm.  If  the  distance  travelled  by  the  ion  is  L  *=  0.8  pm,  which  is 

R 

the  minimum  necessary  distance  needed  to  get  the  ion  completely 

excited  into  the  n»3  states,  then  the  total  energy  loss  is 

<5E  =  NTSLR=e  50  keV.  The  relative  decrease  in  speed  is  given  by 

-  2 

6 v/v  =  6E/2E  ~  4x10  .  These  numbers  show  that  the  constant  velocity 

assumption  is  not  too  bad,  in  the  sense  that  the  error  in  velocity  is 
on  the  order  of  4%  per  Rabi  cycle.  We  should  point  out  these  are  bulk 
estimates,  and  are  therefore  likely  to  be  much  larger  than  the 
corresponding  surface  rates.  As  the  ion  slows  down,  U31  and  become 
detuned.  In  terms  of  the  stopping  power,  the  detuning  from  the  pump 
frequency  is  given  by 


A 
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k  6  v 
n 


vN  t  k  S 
T  R  n 

2  E 


(2.39) 


Estimates  of  ftA0  and  h A  j  are  given  in  Table  I. 

4.  3 

Another  important  assumption  in  the  above  discussion  is  that  Be 

+  3 

ion  remains  Be  inside  the  material  medium,  that  is,  it  keeps  to  its 
original  charge  state.  There  is  some  conflicti’-a  opinion  and  evidence 


on  the  states  of  an  ion  or  atom  while  it  travels  inside  a  solid 
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material  [23-27].  The  projectile  may  lose  its  original  electron  and 
be  completely  stripped,  or  it  may  capture  electrons,  losing  them  only 
while  leaving  the  solid  material  [21].  At  any  rate,  the  charge  states 
of  a  projectile  may  not  be  entirely  determined.  In  view  of  this 
situation,  how  reliable  is  the  above  assumption  and  the  consequent 
analysis?  There  are  three  possible  answers  to  this  question.  First, 
the  charge  state  of  the  ion  is  clearly  dependent  on  the  ion  speed.  At 
low  speeds  the  charge  state  appears  to  remain  unchanged  with 
relatively  high  probability  [24].  One  can,  therefore,  reduce  the 
ionic  speed  by  matching  the  relevant  ion  transition  frequency  to  a 
higher  harmonic  of  the  material  field.  There  is  probably  some  trade 
off  point  in  this  procedure,  since  the  Rabi  frequencies  of  the  higher 
harmonics  of  the  material  fields  are  lower,  leading  to  larger  energy 
losses.  Second,  even  for  swift  ions,  the  ion  spends  some  of  its  time 
inside  the  solid  in  its  original  charge  state  due  to  repeated  capture 
and  loss  of  electrons  of  the  medium.  From  the  cross  sections  for 
electron  capture  or  loss  [14,28,29],  one  can  estimate  the  probability 
of  the  ion  being  in  its  original  charge  state,  and  accordingly  find 
the  expected  values  of  the  desired  quantities  from  this  probability. 
This  picture  of  ions  being  partially  in  their  original  charge  states 
is  supported  by  the  correlations  observed  between  ions  and  a  coterie 
of  electrons  that  are  emitted  secondarily  [30].  Finally  and  third,  in 
the  treatments  of  the  Okorokov  effect,  one  uses  the  picture  of  the 
projectile  being  in  its  entry  charge  state  [7,10,11].  One  then 


calculates  the  probability  of  excitation  of  the  projectile  as  a 
function  of  its  speed.  The  model  apparently  works,  and  there  is  a 
fair  agreement  with  experiment  [12-14]. 

The  cross  sections  for  the  electron  capture  and  loss  are  on  the 
order  of  a c  ~  ~  10' 19  cm2  [14,24,28].  The  target  mobile  electron 

density  for  MnO  is  on  the  order  of  N  ~  10 1 7  cm'3.  It  follows  that 

e 

the  probability  of  the  ion  changing  its  charge  state  i3  roughly 
(o  +  o£)NgLR  ~  10' 6  per  half  Rabi  cycle. 

III.  GAIN  FOR  THE  X-RAY  SIGNAL 


In  order  to  calculate  the  gain  for  the  signal  field,  we  follow 
the  well  known  procedure  used  in  traveling  wave  amplifiers  and  free 
electron  lasers  [31].  We  calculate  the  transverse  charge  current 

i  u)  t 

density  of  the  ionic  electrons  which  oscillates  a3  e  s  ,  and 
substitute  it  into  the  wave  equation.  The  wave  equation  yields  a 
dispersion  relation  which  requires  either  temporal  or  spatial  growth. 
The  latter  yields  the  small  signal  gain  for  the  X-ray  signal. 

The  n-3  n-2  transitions  cause  a  particle  current  density  that 

is  given  by 
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It  follows  from  (2-29)  and  (2-35)  that  the  component  which  oscillates 

i  to  t 

as  e  s  is 
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where 
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Using  (2.24)  and  (2.36),  Dj  becomes 


D1  -41”l‘FlU"))2  * 


(3.4) 


j  i3  the  transverse  current  density  of  just  one  ion.  If  there  are 

s 

N.  ions  per  unit  volume  in  the  ion  beam  moving  with  velocity  v,  then 

f  s) 

the  total  transverse  charge  density  component  is  given  by  ej  ■  N  ^  and 
the  wave  equation  becomes 


“*+  A  c  “  c 

where  A  g  -  yAse  e  is  the  positive  frequency  component  of  the 

signal.  Performing  the  indicated  operations  in  (3.5),  one  obtains  the 
dispersion  relation 


4  0  2 

4Te>rBQ 


N.D 

l 


4 

1 


4  4  4  +  _ 

m  c  k  (ft  -9.  ) 
n 


(n 


1.2 

4*  2 


(SI  -w_-+aj  )t  + 
j2  s' 
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For  an  amplifier,  one  lets  the  frequency  be  real  but  allows  the 
wavevector  to  have  an  imaginary  part  in  order  to  determine  the  growth 
or  the  decay  of  the  signal  in  space.  Thus,  let 


k  = 
s 


u 

s 

c 


is 


(3.7) 


Here  a  positive  g  means  that  the  signal  wave  grows  as  exp(gz)  and  thus 
g  is  the  gain  due  to  the  ions  moving  with  velocity  v.  g  is  generally 
quite  small  compared  to  ~./c;  therefore, 
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(3.8) 


and 
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(3.9) 


In  the  limit  r2  0,  the  terms  bracketed  are  reduced  to  a  pair  of 
delta  functions: 
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In  an  ionic  beam,  ions  will  in  general  have  a  dispersion  of 
velocities  around  an  average  velocity  v0z.  Ions  moving  with  slightly 
different  velocities  will  be  differently  tuned  to  the  material  field 
and  contribute  to  the  overall  gain  with  different  strengths.  Assuming 

+  3 

that  the  beam  of  Be  has  a  narrow  Gaussian  distribution  of  velocities 
in  the  z-direction  described  by 
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(3.11) 


the  net  gain  is  given  by 


(3.12) 


g=  Pcs  f  dv  g(v)  p^v) 


where  we  have  multiplied  the  average  over  the  velocity  profile  with 
pcs  -  the  probability  of  the  Be 3  ion  being  in  the  charge  state  +3 
(roughly  0.5).  If  Tz  <  knAv,  one  can  use  the  simpler  form  of  the  gain 
given  by  (3.10).  The  integral  in  (3.12)  is  then  readily  done  and  one 
f  inds 
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where  vt  is  some  tuning  velocity  given  by 
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(3.14) 


The  expression  (3.13)  for  g  can  be  put  into  a  more  recognizable 
form  in  order  to  compare  it  with  the  usual  small  signal  swept-gain 
[8].  Define  an  "effective"  Einstein  A-coef f icient 
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erf 
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(3.15) 


and  the  Doppler  broadening  half-width 


Aw  =  (k  Av)  2/fn2 


(3.16) 
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(3.13)  can  be  written 
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where  ^nvo*  It  is  seen  that  (3.17)  differs  from  the  usual  gain 

expression  in  two  ways.  First,  instead  of  the  usual  population 

inversion  factor  aN  -  N  ,  -N  .  (3.17)  has  a  density  factor 

excited  ground 

that  refers  to  the  total  density  of  the  ions  modulated  by  the  Rabi 
frequency  and  the  detuning  of  us  from  j32.  In  other  words,  the  entire 
collection  of  ions  participate  in  the  excitation  and  emission 
processes.  This  is  the  origin  of  the  phrase  "coherent  excitation." 
Second,  (3.17)  has  a  Doppler  factor  that  has  complicated  ^-dependence 
and  goes  to  zero  exactly  at  us  «  u>32.  Thus  the  following 
correspondences  exist  when  compared  with  the  usual  gain: 
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and 
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(3.23) 


As  seen  from  Table  I,  g  ^  10  "4  cm'1  for  low  harmonics.  Due  to 

o 

slowing  down  of  ions,  one  can  expect  6 p  to  be  on  the  order  of  An/fiR 
which  is  a  number  generally  larger  than  10.  Because  of  the  factor 
(vq/Av),  it  is  desirable  to  have  as  small  (av/vq)  as  possible.  By 
pushing  to  the  limits  of  the  present  ion  beam  technology  it  is 
possible  to  get  a  high  velocity  resolution  on  the  order  of 
(Av/vq)  10' 5  .  w  is  then  0.04  and  0.15  for  the  zeroth  and  first 
harmonic  frequency  matching,  respectively.  Clearly,  the  maximum  gain 
is  obtained  when  5  is  approximately  given  by 


6  =  -  ~  6  +y  J  6  2+l 

m  2  p  4  p 
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which  makes  the  exponential  factor  unity.  For  5  =  6, 
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Thus,  for  N-  -  10 1 4  cm  i  and  Av/v  -  10  5 ,  g  ^  10  cm  1  which  is  a 

1  o 

significant  gain.  Note  that  the  gain  is  confined  to  a  narrow  tuning 


■  -'■>  -N  . 


range  between  a>  and 

s  J  *- 

expanded  as 


Near  6  the  exponential  factor  can  be 
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exp  % 


(3-26) 


n 

It  follows  that  for  6  -  5  -v  (w2/5  )  which  is  on  the  order  of  10 ' 4  or 

m  P 

less,  g  is  substantially  reduced. 

IV.  CONCLUSION 

With  slight  modifications  the  preceding  analysis  applies  to 

non-hydrogenic  ions  and  neutral  atoms.  We  therefore  conclude  that 

ion-antiferromagnetic  surface  scattering  is  a  promising  and  quite 

general  scheme  for  X-ray  amplification  and  X-ray  lasers. 

Note  that  for  N  -  10’- 4  cm*  3 ,  the  charge  current  density  of  the 
i 

ions  is  JejN.v  Z  t,  6.5x104  amp/cm2  for  v=10?  cm/sec.  Or  in  terms  of 

1  o 

the  power  contained  in  the  ionic  beam,  P.  -  'iN.M  v3  ^  9xl04  MW/cm2. 

ii  ion 

Clearly,  in  an  actual  experiment,  the  particle  beam  would  have  to  be 
pulses  of  a  nearly  neutral  plasma  beam. 

Admittedly  the  beam  densities  required  to  achieve  g  ^  10  cm*1  are 
high  [32].  However,  presently  there  are  sources  for  Li  ions  that  can 
yield  such  high  densities  [33].  Furthermore,  ion  beams  can  be  focused 
on  the  surface  to  obtain  high  densities  [34].  The  experimental 
difficulties  may  lie  with  combining  a  high  velocity  resolution  with 
high  densities. 

Another  possibility  is  to  increase  the  effective  range  of 


interaction  L 


between  the  signal  and  the  excited  ions.  The  actual 


amplification  of  the  signal  amplitude  is  given  by  e  °  .  If  L  ^  is 

sufficiently  large,  much  smaller  g  can  be  tolerated  to  obtain 
significant  signal  amplification.  Lg  ^  can  be  increased  in  a  storage 
ring  configuration,  where  both  the  particle  beam  pulses  and  the  X-ray 
signal  repeatedly  scatter  from  surfaces  that  are  placed  in  a  circular 
fashion  (with  a  large  radius  to  allow  for  grazing  angle  scatterings). 
Between  surface  scatterings,  the  velocity  spread  of  the  particle  beam 
can  be  adjusted  by  means  of  external  fields. 

Finally,  the  preceding  analysis  shows  that  the  antiferromagnetic 
surface  scattering  is  an  extremely  efficient  pump  mechanism  to  create 
an  inverted  population.  which  determines  the  rate  of  production  of 

the  excited  states  is  four  to  five  orders  of  magnitude  larger  the 
radiative  lifetimes  of  the  relevant  excited  states. 
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FIGURE  CAPTION 


Fig.l  The  schematic  of  the  antiferromagnetic  surface  scattering. 
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We  show  that  a  sizable  electric  dipole  moment  can  be  induced  between  the  'if  and  'A,  states  of 
O2  when  oxygen  is  placed  in  an  environment  of  particles  with  spin  -5-. 


I.  INTRODUCTION 

In  diatomic  molecules  formed  by  the  group-VI  elements 
(e.g.,  02,S2,S0)  and  by  a  combination  of  group-V  and 
group-VII  elements  (e.g.,  NF,NC1,PC1),  the  spin  selection 
rule  AS  =0  leads  to  highly  metastable  states  which  can  be 
very  useful  for  high-power  laser  systems.  A  particularly 
well-known  case  is  the  'Ag  state  of  oxygen  molecule.  Ow¬ 
ing  to  spin  and  parity  symmetry,1  its  decay  to  the  ground 
state  is  highly  forbidden  with  spontaneous  decay  rate 
2.6X10-4  sec-1.2  This  fact  is  used  in  oxygen-iodine 
lasers  to  store  the  energy  in  the  ‘Ag  state  of  02,  which  is 
then  collisionally  transferred  to  atomic  iodine  via  the  reac¬ 
tion  02('Ag)  +  I(2/>3/j)  — 02(32-)  +I(2P,/2).  Radiative 
emission  takes  place  on  the  iodine  through  its  decay  back 
to  2P|/2,  which  is  a  magnetic-dipole-allowed  transition 
with  a  spontaneous  transition  rate  7.7  sec-1.3  The  fre¬ 
quencies  of  'Ag-32^  and  2/>1/2-2/>3/2  transitions  match, 
which  makes  it  possible  to  resonantly  transfer  energy  be¬ 
tween  02  and  I.  In  practice  it  is  difficult  to  find  such 
matching  partners  for  other  metastable  systems  which 
might  have  applications  to  high-power  lasers.  For  this 
reason,  as  well  as  to  avoid  losses  involved  in  collisional  en¬ 
ergy  transfers,  there  has  long  been  an  interest  in  obtaining 
radiative  emissions  directly  from  the  excited  metastable 
system  itself.  To  achieve  this  it  is  necessary  to  break  the 
various  symmetries  of  the  metastable  system.  For  in¬ 
stance,  the  parity  symmetry  can  be  broken  by  applying 
static  electric  field  or  nonresonant  radiation,  or  by  means 
of  nonresonant  collisions  with  another  molecular  or  atom¬ 
ic  system.4-6  On  the  other  hand,  breaking  of  the  spin 
symmetry  requires  application  of  either  spatially  varying 
magnetic  field  as  in  a  Stem-Gerlach  device,  or  a  direct 
spin-spin  coupling  to  another  system. 

In  this  paper  we  consider  the  latter  possibility  for  the 
singlet  oxygen,  where  the  electrons  of  02  are  coupled  to  a 
set  of  external  spins  through  their  internal  magnetic  mo¬ 
ments.  This  situation  can  be  realized  by  putting  02(  ) 

in  contact  with  a  magnetic  surface  as  discussed  in  Sec.  III. 
By  means  of  a  simple  perturbation  analysis,  we  show  that 
both  parity  and  spin  symmetries  are  broken  and  that  an 
electric  dipole  moment  is  induced  between  'Ag  and 
which  is  on  the  order  of  10-2  D.  This  implies  a  spontane¬ 


ous  transition  rate  of  about  102  sec-1,  which  is  nearly  6 
orders  of  magnitude  larger  than  the  corresponding  rate  for 
isolated  02. 

II.  INDUCED  DIPOLE  MOMENTS 


Consider  the  situation  illustrated  in  Fig.  I.  To  simplify 
the  problem,  we  take  into  account  just  two  of  the  valence 
electrons  of  02.  These  electrons  are  coupled  to  spin- 7 
particles  at  sites  Rk  (k  — 1,2, . . . ,  TV)  by  means  of  magnet¬ 
ic  dipole-dipole  interactions 

V=z28F-b  2  2  - 7[3 ay{Kk-Xj) 

X=)  j  m  1,2  I  R*  — |  5 


X  ov(Rx-xj) 

-  |R-x;  |2dyo*]  , 


(1) 

where  fiB  is  the  Bohr  magneton,  ok  are  the  Pauli  spin 
operators  of  the  external  spins,  and  g  is  their  gyromagnet- 
ic  ratio.  ct\  and  cr2  are  the  Pauli  spin  operators  for  the  02 
electrons.  In  the  following,  we  assume  that  the  external 
spins  and  02  are  fixed  in  space  and  ignore  all  other  in¬ 
teractions  between  the  two  systems. 

To  first  order,  the  perturbed  states  of  the  combined  sys¬ 
tem  is  given  by 


\a,M)v=\a,M)+  2 

flVfl, 


a\M’) 


|  V  |  a,M) 

Ea  +  —  £,•  —  Wy 


(2) 


FIG.  1.  Diagram  for  the  notation  of  the  coordinates. 
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where  |  a,  M  )  =  |  a  >  ®  |  M),  |  a)  and  £a  designate  the  unperturbed  states  and  energies  of  02,  and  |  M )  and  Wst  the 
unperturbed  states  and  energies  of  the  external  spin  system.  Each  M  refers  to  a  set  of  1 V  parameters, 
|  M )  =  |  . . .  ,sv  )  where  r*  can  take  on  only  two  values  ±1,  corresponding  to  spin  up  and  down  at  the  Ath  site.  It 

follows  from  (2)  that  the  induced  matrix  elements  of  the  dipole  moment  operator  Dop  =  —  e(  X|  +  x2)  are 


y(a,M  |D°P|a',M'>K  =  6MM.2 


(a,M  [  Dop\a",M)(a"M  |  V\a',M)  (  (a,M  |  V  [  a",\f  )(a",M  I  Dop  \a\M) 


Ea-Ea. ■ 


Ea  Eg" 


(3) 


Because  of  the  form  of  the  operator  Dop,  the  sums  over  the  intermediate  states  can  be  performed  in  an  approximate 
way  as  follows.7  Let  us  define 


S(£uii)='£ 

a " 

a" 

such  that 


< a,M  |  V  \a",M)(a",M  I  exp[  —  /(|"|-r1  +  52-F2)]  |  a',M) 


Ea~Ea: 


{a,M  |  exp[  —  i  (|*rT|  +  |*2-r2)]  |  a",M  \V\  a'.M) 


Ea—Ea.. 


v(a,M  1  D op  |  a',M)y=  ~ie 


a  a 
af.  +  af2 


(S+T) 


S,=  S2=o 


(4a) 

(4b) 

(4c) 

(5a) 

(5b) 


Let  us  also  define  the  operators  F  and  G  such  that 

|[F,//o]  +  <£a-£a>£)  |a')=exp[-/(f1-r1  +  |’J-r2)]|a'>  , 

(a  |  j-[G,//0]  +  (£a-£a  )G  =  <a  |  exp[ -«(f,-f1  + f2-r2)]  , 

where  H0  is  the  unperturbed  Hamiltonian  of  02.  One  then  has  S  =  ( a,M  \  VF  |  a',M )  and  T  =  ( a,M  |  GV  \a',M).  To 
evaluate  F  and  G,  we  assume  that  they  depend  only  on  position  variables.  This  approximation  permits  us  to  replace  H0 
in  (5a)  and  (5b)  with  the  kinetic  energy  operator  (pf/2 m  +  p2/2m)  of  the  valence  electrons.  Fand  G  are  then  readily 
evaluated  in  the  plane-wave  representation,  and  S  and  T  become 

<^(k,,k2)Yl(l,2)<A/  I  r(k,-k'1,k2-ki)|M>Ya.(l,2Wa.(k'i  +  f„k'2  +  f2) 


2 

r,.k2.r;.ki 

r=  2 

k,,k2.r;.k  j 


„  ( k '  )2  _  (ic j+fji2 

Ea-Ea-  2  -^—  +  2  -  ,  ■ 

j  =  i,2  ^ m  j  _  i,2  2m 

^(k,-f„k2-f2)^(l,2)(A/  j  Z^(k|  —  k  *|,k2  —  k  2)  \M)Xa\  l,2)d„  (k  ),ic  2) 


Ea-Ea+  2 


j  — 1,2 


2m 


2  ¥ 


(6a) 


(6b) 


7  =  1.2 


2m 


Here  7'"(k1,k2)  is  the  Fourier  transform 


7^(k(,k2)=  J  dxld\1e 


F(x„x2) 


-  _  ,  (oi-k|)(o^-k|) 

3  f 

O  — • 

kj.o 

0  ;  ak  3  _ 

1  k,  | 2 

+6r„^ 


-i  k,R, 


3 


(a2-k2)(dvk2) 


I  k2 1 2 


(7a) 


We  also  separated  the  02  wave  'unctions  into  spatial  and  spin  parts  as  in  i/'„( X|,x*2)Va(  1,2),  and  6a  is  the  Founer 
transform  of  i/>a : 


^(k„k2)=  J  dx,dx2e  *j)i/ia(x,,x2)  . 


(7b) 


It  follows  from  (4c)  that 


MAGNETIC  SURFACE-INDUCED  SPIN-FLIP  TRANSITIONS  IN  . . . 


C.  3 


y(a,M  |D°p|a',Af)>>  =  —  -  '  .  2  ^(k 'i.k  iWfl(k -  k„k  2-k2)(k,  +  k2) 

m(E„  Ea')  r-  r-  r-r- 

k  ,,  k2.  k  k  2 

X[x\n,2){M  |7''(kI,k2)|M)Arad,2)]  . 


Now  let  a  be  the  state  'Ag  with  ML  =  +2,  and  a'  the  state  ilg  with  Ms  =  +  1.  The  02  wave  functions  for  these  states 
are  given  by" 

il>a=irf{x\)irf(xi) ,  (9a) 

*a  =  TT[*+(1)*-(2)_*-(1)*+(2)] ,  (9b) 

^«-  =  -^=-[ir^(x,)ir“(x2)-n-g  (x,)ir/(x2)]  ,  (9c) 

*„•=*+( l)*+(2)  ,  (9d) 

where  ir~(x)  are  molecular  orbitals  for  individual  electrons,  and  X+  are  the  electronic  spinors  for  spin  up  and  down.  A 
straightforward  algebra  then  yields 

D=  y(  x\,Ml  =2’,M  |  D°p  |  }2fMs  =  1  ;M)V 


ilSegp),# 


32  trm(AE)  x 


2<M|oXj|A/)~  /  dx[ff/(x)] 


X-Ril4 


IT.  (X)  , 


where  A £  =£a  -Ea-  is  the  unperturbed  energy  difference 
and  e*  are  the  complex  vectors  e+  —x±iy.  The  matrix 
element  {M  |  \  M )  =Sx-  ^g(x)  can  be  written  in 

terms  of  the  n  =2,  /  =  1,  m  =  ±l  atomic  orbitals  p  +  (x) 
as8 

vf(x)=p±(x  —  a)—  p  +  (x  —  b)  ,  (11) 

where  a  and  b  are  the  positions  of  the  oxygen  atoms.  To 
obtain  an  estimate  for  (10),  we  first  simplify  (11)  by  ap¬ 
proximating  TTg  with 

nj(x)~—rab-Vp±(x)  ,  (12a) 

where  Fa6  =  a  — b  is  the  intemuclear  distance  of  02. 
Second,  we  assume  that  external  spins  are  relatively  far 
from  02  and  use  the  expansion 

1  _ 1  4(xR)  [12(R-x)2  — RV] 

I  x -R  |  R*  Rb  R% 


Using  Eqs.  (12a)  and  (12b),  the  integral  in  Eq.  (10)  can  be 
readily  evaluated.  The  result  is  somewhat  long;  however, 
it  can  be  simplified  by  taking  an  average  over  the  orienta¬ 
tions  of  the  molecular  axis  (i.e.,  over  ?aj, ),  and  also  over 
the  orientations  of  the  external  spins  (i.e.,  over  Rk),  if  N  is 
sufficiently  large  and  the  spins  are  distributed  around  the 
molecule.  One  then  finds 

nSegntfrl,  I 

*  D  I4rrm(i£)r  3  _  +*  ?  ’  "3I 

The  sum  in  Eq.  (13)  can  be  replaced  by  an  expression 
which  involves  the  probabilities  of  finding  the  spin  up 


(P+ )  or  down  (P_ )  at  a  given  site,  and  a  structure  factor 
for  the  external  spin  system, 

<i4) 

X  KX  X  A  A.  r0 

where  r0  is  a  scaling  length  such  that  r0=  10~8  cm.  Z 
may  be  calculated  in  a  similar  manner  to  the  calculation 
of  a  Madelung  constant.  For  an  order  of  magnitude  deter¬ 
mination,  however,  one  may  replace  the  sum  over  A.  by  a 
three-dimensional  volume  integral,  and  Z  is  approximate¬ 
ly  given  by 


z  =  2 


4  7T  r0 

3  r. 


where  r,  is  the  average  distance  between  two  spin  sites 
(hence  r,  is  the  smallest  length  in  the  set  |/?xj).  For 
many  magnetic  materials  rJ/r0~2.9  Setting  g=2, 
A E  —  1,  eV,  and  rab/r0~  1.5,  we  find 

«D»  a(3x  10-3  debye)i(3e_+z)(F+  —P_).  (16) 

In  other  words,  the  induced  moment  is  on  the  order  of 
10“ 2  debye,  and  the  induced  spontaneous  decay  rate  is  on 
the  order  of  102  sec-1.  Furthermore,  if  the  perturbed  02 
is  in  a  laser  field  of  power  p,  then  its  Rabi  frequency  is 

~  1/2 

=  _  10’  sec-’,  (17) 

rt  n  c 

for  p  =  107  W/cm2. 
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III.  CONCLUDING  REMARKS 

It  is  clear  from  the  preceding  discussion  that  if  the 
external  spins  have  no  preferred  direction  and  are  com¬ 
pletely  random,  the  induced  dipole  moment  vanishes 
within  the  approximations  made  above.  In  order  to  obtain 
dipole  moment  between  and  32^~  states,  the  singlet  ox¬ 
ygen  must  interact  with  a  collection  of  particles  which  has 
net  magnetization  over  some  finite  volume. 

The  best  way  to  establish  contact  between  the  singlet 
oxygen  and  a  magnetic  surface  may  be  by  means  of 
grazing-angle  scattering.  We  have  in  mind  the  situation 
in  which  slowly  moving  ( v~  103 —  10s  cm/sec)  singlet  oxy¬ 
gen  molecules  scatter  from  a  ferromagnetic  surface  at 
grazing  angles  (  <  1*),  and  the  stimulated  emission  occurs 
while  molecules  are  in  the  vicinity  of  the  surface.  This 
type  of  scattering  permits  maximum  exposure  of  the 
metastable  molecules  to  the  localized  spins  of  the  magnet¬ 
ic  material,  and  minimizes  the  damage  to  the  surface. 
The  same  surface  can  therefore  be  used  repeatedly.10 

An  important  advantage  of  the  above  scattering  scheme 
is  that  it  cuts  down  on  possible  chemical  reactions  be¬ 
tween  the  singlet  oxygen  and  the  surface  by  restricting  the 
duration  of  the  contact.  In  chemical  reactions  between 
the  singlet  oxygen  and  the  ferromagnetic  materials,  the 
quenching  rate  kg  is  typically  108  mol-1  sec-1.11  Let  us 
consider  a  single  02  molecule  and  assume  that  it  follows  a 
straight  line  trajectory  on  the  surface  with  a  constant  velo¬ 
city  v  for  a  distance  /.  Let  a0  be  the  effective  cross- 
sectional  area  of  this  molecule.  The  number  of  effective 
moles  of  the  material  atoms  that  it  sees  on  this  trajectory 


is  o0lpst/NA,  where  pM  is  the  number  of  atoms  per  unit 
volume  of  the  target  material  and  NA  is  the  Avogadro’s 
number.  Thus  the  overall  transition  rate  is  kga^p^  /NA . 
Since  the  duration  of  the  interaction  is  l/v,  the  probability 
of  quenching  is  given  by 


prQ~ 


^■Q(Jcf2Psi 

Nav 


(18) 


on  the  basis  of  this  simple  model.  Taking  NA  =6x  1023, 
/  =  1  cm,  oq~  10~14  cm2,  pM  —  1022  atom/cm3,  and  u  ~  103 
cm/sec,  one  finds  that  prQ  —  10“  ",  which  rules  out  chem¬ 
ical  quenching. 

There  are,  of  course,  other  types  of  quenching  mecha¬ 
nisms,  such  as  electron  pick-up  or  electron  loss  by  the 
scattering  molecule,  exchange  of  the  molecular  rotational 
energy  with  the  surface,  etc.  Furthermore,  the  surface 
may  partially  be  demagnetized  as  a  result  of  the  scatter¬ 
ing.  However,  we  do  not  expect  these  processes  to  alter 
the  order  of  magnitude  results  presented  above.  We  con¬ 
clude  that  the  scheme  proposed  here  appears  to  be  quite 
promising  to  obtain  an  oxygen  laser.  More  importantly, 
the  scheme  may  be  used  in  obtaining  lasers  from  other 
metastable  species  for  which  collision  partners  are  not 
known  at  the  present  time. 
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V  1  V 


D.  3 


field  are  small.  When  this  is  the  case,  Eqs.  (l)-(3)  may  be 
linearized  to  obtain 


3K2(z,t0,|j0)/3z  =  Kl Cz»  t0 , y0) 


(7) 


*  inoz 

3K1(z,t0,u0)/3z  =  e  Es(z,t) 


(8) 


3Es(z»f)/3z 


=  iD  A. 

i 


r  -i u0z 

I  duo  J(tq  ,  vi0)  e  K2(z,t0,uo) 


(9) 


K2  =  ~(l/2n)  fd 0 o  exp(-i ©0 ) 69  and  f<!  =  -<1/2tt)  d0o  exp(-i0o) (y  -  p0) 
are  respectively  the  density  and  energy  bunching  amplitudes.  We 
expect  the  FEL  to  saturate  when  |K2|  becomes  comparable  to  one.  The 
initial  values  of  K2  and  Ki  are  zero. 

If  there  is  no  inhomogeneous  broadening,  we  set 
JJ(t0,p0)  =  I  (To)  5  (Ho)  •  Then  Eqs.  (7)  —  (9)  may  be  solved  by  Laplace 
transforming  in  z  to  give 


Es(z,t)  =  Es(0,t) 


+ 


~r  FqO)Es(0,t  -  z  72y^c) 


(10) 


where  Fq(x)  is  a  generalized  hyper geometric  function  given  by  the 
series  expansion 


Fo(x) 


E 

n=l 


1 

n !  (2n  -  1) ! 


n 


x 


(11) 


and 
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t-z'/2y^c 

s 


=  i<D| | 22y^cz'  2  dT0  I(x0) 


t-z/2y3c 

s 


Equation  (10)  expresses  the  field  at  position  z  as  the  sum  of  the 
incident  field  at  time  x  and  an  integral  over  the  incident  field  at 
earlier  times.  We  see  explicitly  that  coherence  cannot  develop  in 
one  pass  over  more  than  the  slippage  time  L/2y|c,  where  L  is  the 
wiggler  length.  This  is  true  also  in  the  strong-signal  regime.  In 
other  words,  an  FEL  operating  without  mirrors  by  ASE  cannot  be 
coherent  over  more  laser  periods  than  there  are  wiggler  periods. 

In  the  case  that  both  the  current  and  the  laser  field  are  cw  we 
put  Es(z,t)  =  Es(z)  exp(i2Y^c  Uqx)  in  Eq.  (10)  to  obtain 


fz  1 

Es(z)  =  Es(0)[l  +  dz*  jr 

0 


-iu0z’ 

:'2(z  -  z'))e  ] 


where  C  =  <D|Ai|2I.  Since  the  power  series  (11)  is  very  rapidly 
convergent,  the  field  can  generally  be  calculated  accurately  by 
keeping  only  a  few  terms  in  the  expansion.  Only  the  first  term  is 
needed  to  obtain  the  well-known  antisymmetric  small-gain  formula  with 
peak  gain  .27  CL3  at  y0L  =  2.6. 

An  alternative  approach  to  cw  theory  is  to  drop  the  dependence 
on  x  and  tq  in  Eqs._  (1)  —  (3)  and  (7)  —  (9)  -  In  this  case  we  take  J7(vj0) 
to  be  centered  at  u"0  with  width  U.  lhat  is,  the  detuning  Fo  from 
resonance  is  included  in  J7  rather  than  Es.  If  we  define  the 
normalized  distribution  f ( u o )  by  27  =  If,  then  Eqs.  (7)  —  (9)  become 


dK2(z,y0)/d2  =  , 


*  1U0Z 

dK! (z, u0)/dz  =  kA.  e  E  , 


f  -i^o- 

dEs(z)/dz  =  iD  Ai I  jdy0f(u0)e  K2(z,y0) 
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There  equations  may  be  solved  by  Laplace  transforming  to  give 


a+i°° 

Esiz)  =  2ir  Cs(0)  d8e  Z[B  -  iC  du0f Cuo)/ (B +  ip0) 2]' 1 ,  (17) 

1  ^  J  a-i= 

where  the  contour  is  to  be  taken  to  the  right  of  all  singularities. 
The  simplest  case  is  where  f(Po)  is  a  Lorentzian, 


Then  the  integral  over  yg  in  Eg-  (17)  gives  (8  •+  iy 0  +  U)2.  We  may 
close  the  contour  integral  over  8  around  the  three  poles  at  the  roots 
6^  of  the  equation 


8(8  +  iy0  +  U)2  =  iC 


(19) 


and  evaluate  the  integral  to  express  Es(z)  as  a  linear  combination  of 
the  three  exponentials  exp(8jz) .  In  the  case  U  =  0  it  can  be  shown 
that  this  solution  is  equivalent  to  that  in  Eq.  (13) .  In  general  at 
most  one  of  the  three  roots  has  a  positive  real  part,  and  the 
corresponding  mode  will  dominate  near  the  end  of  a  sufficiently  long 
wiggler  (unless  saturation  sets  in  first:)  .  The  detuning  giving 
maximum  growth  of  the  unstable  mode  is  uo  =  U/3  -.  The  gain  is 
greatly  reduced  if  U  £  Cl/3.  In  the  case  that  U  «  C1/3and  Ug  =  U/3^2 
the  three  modes  contribute  about  equally  at  the  entrance  of  the 
wiggler  and  the  gain  in  the  limit  of  very  large  gain  approaches 

1  U,  1/.  4 

G  =  9  exp[(3"c  -  3  U)L]  .  (20) 

Graphs  showing  the  dependence  of  2Re(8)  on  yg  and  U  may  be  found  in 
Reference  4.  The  condition  that  the  FEL  saturates,,  when  K ?  =  1  can  be 
used  to  infer  a  saturation  power  Psat  =  (l2A^/4k  )73(c/eoE)Mi. 

SEMICLASS ICAL  FEL  GAIN 

When  electrons  emit  a  laser  photon,  they  undergo  a  recoil 
6E  If  we  let  the  recoil  in  energy  detuning  units  be  2q,  we 

see  from  Eq.  (4)  that 


D. 


q  =  hk2/2Mcy|  . 


(21) 


If  2q  is  comparable  to  or  larger  than  the  homogeneous  broadening  tt/L, 
then  the  classical  gain  formula  becomes  incorrect  and  seriously 
overestimates  the  true  gain.  The  quantum  recoil  can  be  important  for 
an  electromagnetic  wiggler  at  infrared  wavelengths.  In  this  case  y 
is  typically  v.  50  and  i  -  1  is  small.  Note  that  L  for  an  } 

electromagnetic  wiggler  is  only  half  the  wiggler  pulse  length  (L  =  — 
cTp) ,  since  the  wiggler  pulse  and  the  electrons  pass  through  each 
other.  The  product  2qL  may  be  written  as 


2qL  =  582" 


T  (nsec) 

_JL 


[Xi(ym)]%[Xs(A)J,/i 


(22) 


As  an  example,  for  a  C02  laser  pulse  1  nsec  long  used  to  generate 
radiation  of  =  lOOA^Eq.  (22)  gives  1.6". 

It  has  been  shown  by  Renieri6  that  the  condition  2qL  >  ir  implies 
that  the  quantum-mechanical  spreading  of  the  wave  packet  describing 
the  electron  is  larger  than  the  laser  wavelength.  Ths  is  another 
indication  that  a  quantum-mechanical  treatment  of  the  electrons  is 
required . 

We  present  here  a  semiclass ical  theory  of  the  FEL  in  which  the 
electrons  are  described  quantum  mechanically,  while  the  field  is 
treated  classically.  This  approach  accounts  satisfactorily  for  the 
quantum  recoil  but  (as  in  Lamb's  sariiclassical  theory  of  the  laser6) 
i  yes  not  include  the  noise  needed  to  produce  laser  start-up  when  no 
laser  field  is  initially  present.  In  the  laboratory  frame  the 
semiclassical  theory  can  be  developed  by  coupling  the  Klein-Gordon 
equation  for  the  electrons  with  the  Maxwell  equation  for  the  field. 

An  equivalent  approach  is  obtained  by  recalling  that  0  and  y  are 
proportional  to  the  position  and  momentum  of  an  electron  in  a  frame 
moving  at  the  speed  of  the  pondercmotive  potential.  In  quantum 
theory  these  become  operators  0  and  y  whose  commutator  is 


[0 1  y ]  =  2iq  .  (23) 

Restricting  our  attention  to  cw  operation,  Eqs.  (1)  and  (2)  become 
tiie  operator  equations 


d0/dr  =  y, 


(24) 
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dy/dz  =  -<[A^  Es(z)el9  +  c.c.] 


(25) 


We  consider  here  only  the  case  of  a  uniform  wiggler,  so  A-  is 
constant.  For  a  wiggler  field  which  is  a  high-power  Gaussian  beam, 
the  wiggler  is  effectively  tapered  by  the  slowing  of  the  electrons  in 
the  vicinity  of  the  beam  waist.4  This  is  an  important  source  of 
additional  homogeneous  broadening  if  the  pump  laser  power  exceeds 

neom2c5/2e2  =  1GW. 

Equations  (24)  and  (25)  may  be  interpreted  as  Heisenberg 
equations  of  the  "Hamiltonian" 


H  =  y2  ~  i<A^  Ese10  +  i*Aj.Es  e  (26) 

with  the  following  rule  for  the  evolution  of  operators: 

3A/3Z  =  (i/2q)[H,A]  .  (27) 

Similarly  the  Schrodinger  equation  for  state  vectors  becomes 


H|^>  =  2iqd  J  i^>/dz  .  (28) 

This  formulation  is  equivalent  to  nonrelativistic  quantum  mechanics 
in  the  moving  frame.  The  semiclassical  generalization  of  Eq.  (3)  is 
obtained  by  replacing  the  averages  on  the  right  side  by  the  quantum 
expectation  value,  yielding 

dEs(z)/dz  =  DA.I  tr(pe"10)  ,  (29) 

where  p  is  the  density  matrix.  If  we  assume  that  the  initial  density 
matrix  is  diagonal  in  the  u  representation. 


P  (u  0  •  VJQ  )  =  fOo^Oo  -  Uq)  •  (30) 

we  can  use  first-order  perturbation  theory  to  obtain  a  quantum 
generalization  of  Eqs.  (14) - (16) .  It  turns  out  that  only  Eq.  (15)  is 
altered.  It  becomes 
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i 
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*  = 

dKi (z,u0)/dz  =  <Aa  e  Es  '  q2K2  •  (31) 

If  f(uo)  is  Lorentzian,  the  cubic  dispersion  relation  (19)  is 
generalized  to  become 


r- 


B  [  (B  +  iu0  +  U)2  +  q2]  =  ic.  (32) 

We  show  in  Reference  4  that  the  quantum  recoil  acts  to  reduce  2Re(B) 
for  the  unstable  mode  below  its  classical  value,  as  well  as  to  reduce 
the  range  in  m~q  over  which  there  is  an  instability. 

For  the  case  of  small  gain  per  pass,  the  gain  can  be  written 
quite  generally  as 


G  =  C  ~— 

2q 


dp0[f(uo  +  q)  -  f(uo  -  q)l 


u0/2 


(33) 


k 


If  2q  is  very  small  compared  to  either  the  homogeneous  broadening  tt/L 
or  the  inhomogeneous  broadening  U,  then  G  reduces  to  its  classical 
value.  The  condition  U  <  2q  is  rather  restrictive;  it  is  equivalent 
to  the  condition  that  the  energy  spread  of  the  electrons  be  less  than 
huis.  For  100A  x  rays  this  is  less  than  1  keV.  Even  if  low  energy 
electrons  are  used  (as  would  be  the  case  with  an  optical  wiggler) , 
the  relative  energy  spread  needed  to  produce  quantum  modifications  in 
the  gain  would  still  be  less  then  .01%.  However,  even  in  the  regime 
U  >  2q,  if  2q  >  tt/L,  quantum  effects  should  be  manifested,  for 
example,  in  the  statistics  of  the  emitted  photons  or  in  the  electron 
energy  distribution  at  saturation.  When  2q  >  tt/l,  one  expects  the 
FEL  to  saturate  when  one  photon  has  been  emitted  per  electron. 

The  validity  of  first-order  perturbation  theory  depends  on  the 
smallness  of  the  product  of  the  interaction  energy  with  the 
interaction  time  divided  .by  h.  In  terms  of  the  Hamiltonian  (26)  one 
requires  K]AiEs|L  <  2q.  Substituting  for  the  values  of  <  and  q,  this 
becomes  (e*- 1  A^Eg  |/2Mc^Yshks)  L  <  1.  For  an  x-ray  FEL  (in  contrast  to, 
say, an  infrared  FEL)  this  condition  will  hold  for  relatively  large 
laser  intensities  (Es  of  the  order  of  106V/m  with  an  optical 
wiggler) ,  making  first-order  perturbation  theory  appropriate  to  deal 
with  at  least  the  small-signal  regime.  In  the  classical  regime 
(2qL  <<  tt)  it  turns  out  that  the  gain  predicted  by  perturbation 
theory  is  correct  even  for  conditions  where  perturbation  theory  would 
be  expected  to  break  down,  but  it  remains  to  be  seen  whether  this  is 
true  of  the  quantum  regime. 


Further  discussion  of  the  quantum  regime  for  an  FEL  with  an 
infrared  laser  pulse  as  the  wiggler  may  be  found  in  Reference  4, 
where  numerical  examples  are  given  for  which  the  gain  per  pass  may  be 
as  high  as  40%  at  5.7A. 

CONSTRAINTS  ON  A  LONG-WIGGLER  X-RAY  FEL 

One  can  anticipate  that  the  extension  of  FEL  operation  to  x-ray 
wavelengths  by  using  high-energy  electrons  and  a  long  conventional 
wiggler  will  be  a  gradual  and  evolutionary  one,  in  which  one  builds 
oh  the  experience  gained  by  previous  experiments  at  longer 
wavelengths.  Electron  beam  quality  appears  to  be  the  greatest 
problem,  although  there  are  also  difficult  problems  in  resonator 
optics  and  wiggler  construction.  For  soft  x-rays  storage  rings 
provide  an  acceptable  beam  quality,  and  experiments  using  a  storage 
ring  to  drive  a  100A  FEL  are  being  planned.8  Below  about  60A  the 
storage-ring  approach  becomes  untenable  because  of  the  quantum  energy 
spread  in  the  electron  beam  induced  by  synchrotron  radiation  in  the 
ring.4’9  In  principle  one  could  get  a  better  quality  beam  at  high 
energy  by  using  a  linear  accelerator,  but  suitable  accelerators  have 
not  yet  been  built.  It  appears  that  rf  linacs  give  too  much  energy 
spread,  but  induction  linacs  may  offer  a  possibility. 10 

We  now  list  several  constraints  on  the  operation  of  a  high-gain 
FEL,  referring  to  a  numerical  example  from  Reference  4,  the 
parameters  of  which  are  surmiarized  in  Table  1.  Many  of  these 
constraints  may  be  expressed  as  limitations  on  the  length  L  for  which 
a  cold-beam  classical  one-dimensional  gain  calculation  is  correct. 

Table  I  Parameters  for  a  5  A  X-Ray  FEL 


Quantity 


Symbol 


Value 


Wiggler  wavelength  Xq 

Magnetic  field  B 

Mass  shift  A 

Laser  wavelength  Xs 

Electron  energy  E 

Current  I 

Laser  mode  area  l 

Wiggler  length  L 

Gain  G 

Saturation  power  Psat 

Slippage  time  L/^y^  c 

Characteristic  gain  length  3'/3C  ^ 

Filling  factor 

Electrons  per  laser  wavelength  IXs/ec 


3.2  cm  . 

.24  T 
1.512 
5  A 

3.554  Gev 
10  A 
.16  mm2 
276  m 
1000 

1.7  HW_  ' 

1.44  x  10  14sec 
30.29  m 
1 

104 


A.  Reciprocal  quantum  width:  */2q  =  ^McY^/l’k2  =  11.5  km. 
Since  this  is  »  L,  the  device  is  completely  classical. 
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B.  Diffractive  spreading  of  laser  beam:  Rayleigh 

range  =  E/Aq  =  320  m.  This  is  of  the  same  order  as  L,  so  diffraction 
is  a  significant  effect. 

C.  Reciprocal  energy  spread:  1/U  =  ( Xq/4 tt)  (E/6E)  .  If  we  take 
the  energy  spread  6E  to  be  lOkeV  we  obtain  1/U  =  900  m,  which  is 
larger  than  L.  Equation  (20)  indicates  that  the  gain  is  reduced  from 
1000  to  664.  The  difficulty  is  whether  one  can  obtain  such  a  small 
energy  spread  at  this  high  energy.  This  is  presently  unclear. 

D.  Wiggle  amplitude:  One  must  check  that  the  wiggle  amplitude 
does  not  exceed  the  assumed  cross-sectional  dimension  of  the  beam. 

In  the  present  case  the  wiggle  amplitude  is  eBA2c/4Tt2E  =  .5  pm,  so 
there  is  no  problem.  ^ 

E.  The  characteristic  distance  over  which  the  plectron  beam 
spreads  due  to  Coulomb  repulsion11  is  (4q)mc3Y3£/Ie)  z  ~  5.4  km.  This 
is  clearly  no  problem. 

F.  The  characteristic  distance  over  which  the  electron  beam 
spreads  due  to  emittance  e  is  I/ne.  If  we  assume  that  the  normalized 
emittance  en  ~  e  is  the  same  as  for  the  Stanford  superconducting 
linac,  then  e  =  .06  mm  mrad  (43MeV/3.554  GeV)  .  Then  E/tte  =  70  m. 

This  is  small  enough  to  be  a  problem.  However,  it  is  a  problem  which 
should  be  correctable  by  placing  quadrupole  focusing  elements  along 
the  wiggler . 

G.  The  condition  that  electrons  with  emittance  stay  in  phase 
with  the  ponderomotive  potential  imposes  a  constraint  on  the  length 
of  the  wiggler.  The  phase  angle  change  for  electrons  of  velocity 
component  v_^  is 


e/V  | 2/m2c2) ]  , 


(34) 


where  P  is  the  transverse  canonical  momentum. 

Since  kq  -  (kc/2y2) (1  +  e2A|/m2c2)  =  k  -  ks/2y2=  0,  only  the 
terms  P  •  A^  and  P2  survive  in  the  integrand.  However,  P  •  Aj^  is 
oscillatory  and  averages  to  zero.  Therefore 


|  60 |  =  (k  / 2y 2 } (P  2/ m2c2) L  =  (k  P2/m2c2)L  . 
s  q 


(35) 


We  may  relate  P  to  the  emittance  by 
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P2  =  (e2/E)p2  =  (e2/E)m2c2Y2  =  m2c2e2/E  , 
z  n 

so  that  if  we  require  j<^0|  <  2^,  we  get  L<EA  /e2  =  201  m.  This  is 
small  enough  to  be  significant. 

H.  Transverse  dependence  of  the  wiggler  field  leads  to 
dephasing  of  electrons  which  sample  the  field  at  different  transverse 
positions  r.  Let  us  assume  that  the  wiggler  field  varies 
quadratically  off  axis,  such  that  <5A  /A  =  r2/A 2.  The  actual 
electron  motion  is  rather  complicated  (betatronM oscillations) ,  but  we 
can  get  a  rough  estimate  by  letting  A^  vary  with  r  in  Eq.  (34)  (now 
setting  P  =  0)  .  This  gives  |69|  =  kqL5 (e2 j A? | 2/m2c2) .  If  we  assume 
that  e2|Ai|2/m2c2  is  of  the  order  of  unity,  which  is  typically  the 
case  for  magnetostatic  wigglers  used  in  FEL's,  and  let  the  variation 
in  r2  be  E,  then  |  60  J  =  kqLE/A2.  The  condition  1 6 0 1  <  2ir  now  becomes 
L  <  A3/E  =  205  m.  This  is  small  enough  to  be  a  problem.  Note  that 
the  length  constraint  here  is  proportional  to  1/E,  in  contrast  to  the 
constraints  in  items  B,  E,  F  and  G,  which  are  proportional  to  E. 

I.  To  operate  without  a  resonator  by  ASE,  there  must  be 
sufficient  noise  to  get  the  lasing  started.  For  conventional  lasers 
there  is  a  characteristic  delay  time  after  which  the  lasing  starts. 
For  the  FEL  there  is  an  analogous  delay  distance.  According  to 
Bonifacio,  Narducci  and  Pellegrini,12  the  delay  distance  z  is  given 
by  the  product  of  the  characteristic  gainlength  3 *'^C  V*  andu 

(In | N |  +1),  where  N  is  the  number  of  electrons  per  laser  wavelength. 
This  criterion  gives  zD  =  101m,  which  means  that  lasing  should  occur 
without  a  resonator,  at  least  on  most  shots. 

A  disadvantage  of  operating  without  a  resonator  is  that  the 
temporal  coherence  is  limited  to  the  slippage  time,  which  is  only 
1.44  x  10"14sec  for  the  present  example.  Perhaps  an  order  of 
magnitude  increase  in  temporal  coherence  can  be  obtained  by  using  a 
Bragg  crystal  resonator,  simply  because  of  the  narrow  wavelength  band 
6 A  over  which  the  resonator  is  a  good  reflector.13  This  will  extend 
the  optical  wave  packets  to  a  minimum  length  A|/c6A  after 
reflection.  Nevertheless,  it  is  unlikely  that  an  x-ray  FEL  can 
approach  the  coherence  obtainable  by  free-electron  harmonic 
up-conversion  of  an  already  coherent  laser  signal. 

In  conclusion,  it  appears  that  extension  of  FEL  operation  down 
to  the  5  A  range  will  be  difficult,  but  maybe  not  impossible.  It 
will  depend  on  experience  gained  by  doing  high-gain  FEL  experiments 
at  longer  wavelengths  and  on  developments  in  accelerator  technology 
motivated  by  FEL  applications.  An  FEL  with  a  long  conventional 
wiggler  is  capable  of  higher  gain  than  one  using  an  infrared 
electromagnetic  wiggler.  On  the  other  hand,  the  latter  type  of 
device  is  more  compact,  uses  low-energy  electrons,  and  presents  an 
opportunity  to  explore  the  quantum  regime  of  the  FEL. 
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Abstract 

We  present  the  nonrelativistic  theory  of  nuclear  beta-decay  of  arbitrary 
forbiddenness  in  the  presence  of  an  external  plane  wave  field.  Emphasis  is  on 
the  question  of  whether  the  field  is  able  to  break  the  common  selection  rules, 
so  that  forbidden  decays  could  be  significantly  enhanced.  It  turns  out  that 
while  this  is  true  in  principle  the  required  field  strength  is  of  the  order  of 
the  critical  one  which  is  far  beyond  experimental  possibilities.  Recent  claims 
to  the  contrary  are  contested.  In  particular,  in  the  case  of  the  first  order 
forbidden  decay  of  '  Sr,  we  obtain  an  enhancement  which  is  smaller  than  a 
previously  published  estimate  by  twenty- four  orders  of  magnitude. 


I .  Introduction 


Uuclear  decay  rates  are  usually  considered  to  be,  on  earth  and  under  even 
extreme  laboratory  conditions,  almost  on  the  same  footing  with  fundamental 
constants,  rcr  example,  they  are  commonly  supposed  to  be  essentially 
independent  of  changes  in  their  environment  as  specified  by  temperature, 
external  electromagnetic  fields,  etc. 

Recently,  however,  it  has  been  claimed  [1,2]  that  total  decay  rates  of 

forbidden  nuclear  beta-decay  can  be  significantly  enhanced  by  an  intense  plane 

wave  field.  The  idea  was  that  external  fields  be  able  to  change  the  selection 

rules  of  the  decay  or,  at  least,  to  generate  contributions  to  the  decay  rate 

which  are  of  the  same  type  and  order  of  magnitude  as  those  which  usually  allow 

for  forbidden  decays.  Whereas  it  is  beyond  any  doubts  that  this  can  be  achieved 

with  fields  of  the  order  of  the  critical  field  strength  Ecr)li  =  (m^/et) 
lb 

—  1.3x10  Volt/cm,  the  claim  is  that  experimentally  accessible  fields  are 
already  sufficient  in  order  to  render  the  effect  observable,  e.g.,  focused 
intense  laser  fields  with  a  field  strength  of  about  10  10  Volt/cm  or  radio  fields 
of  a  transmission  line  with  even  much  smaller  field  strengths.  The  original 
estimates  [2]  suggesting  enhancements  by  many  orders  of  magnitude  for  highly 
forbidden  decays  were  shown  to  be  erroneous  [3j.  However,  in  Refs.  4  and  5, 
significant  enhancements  are  still  predicted.  The  present  authors  contested 
these  statements  on  the  basis  of  general  arguments  in  Ref.  3.  However,  neither 
in  Ref.  3,  nor  in  Refs.  4  and  5  were  quantitative  estimates  of  the  total  decay 
rate  given. 

In  this  paper,  we  derive  and  evaluate  explicitly  the  total  rate  for  a  decay 
of  arbitrary  forbiddenness  in  the  presence  of  an  intense  plane  wave  field.  Our 
theoretical  model  will  be  essentially  that  of  Refs.  1  and  2.  The  electron  will 
be  described  by  the  Volkov  solution  for  a  charged  particle  in  the  presence  of  a 


plane  wave  field.  The  nuclear  model  will  be  the  shell  model  with  an  inert 
0*-core,  which  is  affected  neither  by  the  decay  nor  by  the  external  field,  and 
one  or  more  valence  nucleon  in  an  angular  momentum  coupled  state,  so  that  the 
nuclear  wave  functions  are  specified  by  an  effective  one-particle  Hamiltonian. 
However,  as  will  be  pointed  out  later,  we  disagree  with  Refs.  1  and  2  on  the 
interpretation  of  the  nuclear  wave  functions  in  the  external  field.  We  shall 
achieve  great  formal  simplification  by  restricting  ourselves  to  a 
nonrelativistic  treatment  in  the  spirit  of  Fermi's  fundamental  paper  [6]. 

Hence,  our  results  are  not  immediately  applicable  to  real  decays.  However,  the 
purpose  of  this  note  is  to  give  simple  and  explicit  results  for  an  albeit 
idealized  situation,  which,  if  positive,  would  encourage  a  more  realistic 
treatment.  Moreover,  we  have  initiated  an  analogous  relativistic  treatment  [3], 
the  conclusion  being  that  the  field  dependence  is  much  the  same  in  the 
relativistic  and  nonrelativistic  case. 

It  should  be  mentioned,  that  we  shall  entirely  disregard  effects  which  are 
due  to  the  atomic  environment.  It  is  well  known  that  external  parameters  such 
as  temperature  do  affect  observed  nuclear  decay  rates  via  their  impact  on  the 
electron's  wave  function,  in  particular,  if  it  is  bound,  like  in  electron 
capture,  bound  state  beta-decay  or  internal  conversion.  These  effects  have  been 
estimated  and  been  found  to  be  of  relative  order  of  10  ,  at  best,  under 

favorable  conditions;  for  a  review  see  [7],  see  also  [8].  In  contrast,  the 
effects  predicted  in  [1,2]  are  much  larger. 

In  Section  II  we  shall  set  up  the  general  framework.  In  contrast  to  Kefs. 
1,  2  and  4,  5  we  shall  directly  calculate  the  total  decay  rate.  This  means,  we 
first  square  the  S-matrix  element  for  the  decay,  then  integrate  over  the  final 
phase  space,  and  do  the  integrations  over  time  last.  This  procedure  which  has 
previously  been  found  to  be  advantageous  in  external  field  problems  (see,  e.g.. 
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Ref.  9)  leaves  the  dependence  on  the  external  field  transparent  throughout  the 
calculation  and  bypasses  the  cumbersome  and  error  prone  integration  over  the 
differential  decay  rate.  In  Section  III  we  proceed  to  an  explicit  evaluation  of 
the  total  decay  rate  for  a  monochromatic  plane  wave  field  with  circular 
polarization  studying  the  possible  enhancement  of  decays  of  arbitrary 
forbiddenness . 

Z 

In  any  case,  enhancements  will  be  seen  to  be  of  the  order  of  (e/E^.^  ) 
with  E  being  the  applied  electric  field  strength.  This  renders  the  effects 
entirely  unobservable,  given  the  limited  field  strength  that  can  for  a  limited 
time  be  obtained  under  laboratory  conditions.  In  particular,  under  the 

conditions  envisioned  in  Ref.  2  the  enhancement  will  be  found  to  be  of  the  order 

— 2V- 

of  10  .  In  concluding,  we  shall  discuss  the  possible  reason  for  the  huge 

discrepancy  between  Refs.  1,  2  and  our  results. 

II.  BASIC  FORMALISM 

In  order  to  make  this  paper  self-contained  and  to  introduce  our  notation  we 
will  here  outline  the  basic  formalism  which  has  been  more  extensively  presented 
in  a  previous  paper  [3j.  We  shall  exclusively  employ  the  Coulomb  gauge  and 
adopt  the  long  wave  length  approximation  for  the  electromagnetic  field  so  that 
A  =  A(t),  in  accordance  with  our  nonrelativistic  treatment.  We  use  units  such 
that  ti  -  c  =  1,  and  the  electron  charge  is  e  =  —  I e I  -  To  lowest  order  in  the 
weak  interaction  gV  and  to  all  orders  in  the  electromagnetic  interaction  the 
S-matrix  element  for  nuclear  beta-decay  is 
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is  the  wave  function  of  the  neutrino  with  momentum  q,  and 

■  e*/><r< J  t  &  -£t 
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with  E  =  p  /2m  is  the  nonrelativistic  Volkov  solution  describing  the  electron 

.r^ 

with  momentum  p.  T  and  T  are  the  initial  and  final  nuclear  wave 

r**  """ 

functions,  which  are  solutions  of  a  one  particle  Schrodinger  equation 


/T 


(2.4) 


with  a  self-ccnsistent  nuclear  potential  V(r)  and  reduced  charges  and  masses  [lj 

.w  J 

and  ,  respectively.  Following  Ref.  1,  we  shall  replace  the  nuclear 
wave  functions  by 
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with  (p  a  solution  of  the  Schrodinger  equation  with  the  Hamiltonian 


-Hc  -  A  ■!.  V ly-)  .  (2-6) 

wv/ 

The  wave  function  (2.5)  would  be  exact  if  were  derived  from  the  Hamiltonian 


H.  -  erE(t)  with  E  =  ~^A/c)t,  rather  than  from  Eq.  (2.6).  It 

*  A/  %/ 


is  important  to 


note  that  the  wavefunction  (2.5),  as  it  stands,  does  not  include,  contrary  to 
appearance , any  interaction  with  the  external  field,  but  is  rather  the  correct 


noninteracting  wave  function  in  the  Coulomb  gauge.  It  is  here  that  our 

interpretation  deviates  from  Ref.  1,  where  it  is  assumed  that  the  ^ ^ 

—  o 

represent  a  fair  approximation  to  the  interacting  nuclear  wave  functions.  V.’e 
have  previously  discussed  this  point  at  length  [3].  Here  we  shall  be  content 
with  pointing  out  that  the  operator  p  =  -  ivin  Eq.  (2.6)  represents  the 

A/ 

canonical  momentum  p  =  Mv  +  eA,  which  does  not,  since  we  use  the  Coulomb  gauge 

~  ~ — 

with  A  f  0,  agree  with  the  mechanical  momentum.  Hence  an  eigenfunction  of  H^ 
does  not  appropriately  describe  the  state  in  the  absence  of  the  field.  This  is 
evident  from  the  expectation  value  of  the  velocity  in  an  eigenstate  of  H^, 

<h)C  //»>  -  <§,  /£,>, 

which  does  depend  on  the  field,  whereas  the  expectation  value  in  a  state  (2.5), 


1  '  t  >  =  <§.  J-i?  /  &  > 


does  not,  as  it  should  be  for  a  noninteracting  state.  Hence  one  has  to  be  aware 


that  by  using  Eq.  (2.5)  one  entirely  ignores  the  interaction  between  the  nucleus 
and  the  external  field.* 


The  total  decay  rate  per  'unit  time  is  obtained  from 
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The  standard  procedure  is  now  to  first  evaluate  the  differential  decay  rate 
IS  .  I  and  then  to  do  the  integration  over  phase  space.  This  would  require  to 
specify  the  external  field  A(t)  from  the  outset  so  that  the  integration  over 
time  in  S,,  can  be  carried  out.  The  final  integration  over  phase  space  is 

■A 

formidable  [1,4,5]  and  has  not  been  done  yet  in  Refs.  4  and  5.  Since  we  are  not 
interested  in  the  differential  decay  rate,  we  shall  here  reverse  the  o».der  of 
integrations:  we  shall  first  integrate  over  phase  space  expressing  the  total 
decay  rate  in  terms  of  the  electron's  and  neutrino's  Green  functions,  and  do  the 
integrations  over  time  afterwards.  The  integration  over  r  in  S.,  need  not 

~  f-L 

concern  us  here;  it  will  be  absorbed  into  a  nuclear  matrix  element  which  is  left 
untouched.  Hence  we  write 


•k 


In  view  of  the  fact  that  the  frequency  and  field  strength  of  the  external 
field  are  very  small  on  a  nuclear  scale  one  would  expect  this  to  be  at  least 
a  fair  approximation.  Unfortunately,  no  exactly  solvable  model  seems  to  be 
available  in  order  to  test  this  assumption.  The  only  exception  is  the 
harmonic  oscillator  which  can  be  solved  exactly  in  the  presence  of  a  periodic 
external  field.  This  solution  corroborates  the  above  assumption.  However, 
the  harmonic  oscillator  is  not  a  very  realistic  example,  given  the 
equidistance  of  its  levels  and  the  lack  of  a  continuum.  So  the  question  of 
whether  the  field-nucleus  interaction  might  produce  noticeable  effects, 


remains  open  to  some  extent,  in  particular,  if  the  nucleus  has  a  very  low 
lying  excited  state  or  if  he  nuclear  energy  release  in  the  beta  decay  is  very 
small.  In  any  event,  to  improve  on  the  approximation  (2.5)  will  be  very 
hard. 
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In  the  preceding  equation  we  introduced  -  E^,  -  m  as  the  nuclear  energy 

release  minus  the  electron  mass  and  have  taken  the  charge  difference  [lj  eT.  - 
to  be  Bi  as  required  by  charge  conservation.  The  Green  functions  are 

-  /4  tyf&Ji? 
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where 
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In  Eq.  (2.12)  we  introduced  the  free  electron  Green  function 
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The  Green  functions  allow  for  explicit  representations  as  exhibited  in  the  last 
lines  of  tqs.  (2.10)  and  (2.14);  however,  the  preceding  integral  representations 


with  respect  to  energy  will  turn  out  to  be  more  convenient. 
Noting  that 
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we  can  now  rewrite  Eq.  (2.9)  as 
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This  representation  has  some  pleasant  features.  The  entire  field  dependence  is 

contained  in  the  functions^  (t’,t)  andij7£(t’ ,t)  so  that  by  dropping  them  we 

recover  the  free  decay  rate  (for  an  arbitrary  degree  of  forbiddenness).  The  • 

quantity  ~$l  is  more  than  just  a  convenient  abbreviation.  It  is  the  only  term 

which  depends  both  on  the  field  and  on  r  and  r*.  Hence  it  is  the  prime  .  .  \ 

candidate  for  a  possible  field  induced  removal  of  fcrbiddenness .  However,  Eq.  ~  • 

(2.15)  shews  remarkable  cancellations:  all  the  individual  terms  on  the  left 

hand  side  are  of  the  order  of  elAlR  with  R  the  nuclear  radius.  For  the 

r>*  U  0 

experimental  conditions  envisioned  in  the  introduction  this  can  readily  be  of  • 

order  unity  or  even  larger.  In  contras  ,  the  right  hand  side  involves  the  .. 
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elec 


trie  field  E  rather  than  the  vector  potential,  and  we  have 
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If  we  assoc, e  that  the  integrals  over  tine  receive  significant  contributions  cnly 
for  Itl  =  |t'-  t|  ^  as  suggested  by  the  uncertainty  relation,  we  have  in  Eq. 
( 2 . 18 )  Uj  T  k-'/Ej  <<  1  so  that  the  impact  of  the  function  Iff,  ( r‘  ,t)  should  be 
marginal.  Even  if  we  do  not  invoke  the  uncertainty  relation,  we  notice  from. 

Eqs.  (2.10)  and  (2.14)  that  the  integrand  of  Eq.  (2.17)  is  proportional  to 
r  <  Notably,  this  rapid  drop-off  is  independent  of  the  external  field.  This 
fact  makes  it  very  unlikely  that  the  integral  (2.17)  will  receive  any 
significant  contributions  from  larger  values  of  T.  With  regard  to  the  second 
field-dependent  function  in  Eq.  (2.17)  we  note  its  power  series  expansion 
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This  starts  with  the  cubic  power  of  Tr,  so  that  it  is  not  likely  to  play  any 
significant  role,  either. 

The  preceding  arguments  suggest  already  that  any  field  induced  effects 
should  be  fairly  small.  Yet,  hand  waving  estimates  of  the  double  integral 
(2.17)  might  seem  to  be  unreliable.  Because  of  that,  and  in  order  to  get 
explicit  numbers,  we  shall  turn  to  an  explicit  evaluation  in  the  next  Section. 
III.  EXPLICIT  CALCULATT I CN  OF  THE  TOTAL  DECAY  RATE  FOR  CIRCULAR  POLARIZATION 
We  shall  now  evaluate  Eq.  (2.17),  explicitly  for  a  monochromatic  circularly 
polarized  field 
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With  slightly  inore  algebra,  linear  polarization  could  be  treated  along  the  sa; 
lines.  The  field  dependent  functions  (2.13)  and  (2.15)  are  new 
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ffote  tha.tTlZ^S  ,t)  does  only  depend  on  and  the  integrand  of  Eq.  (2.17)  only 
depends  on  T  =  (t  +  t‘)/2  via  the  quantity  (t  ,t).  The  integral  over  ~c'  can 


now  be  done: 
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where  r  =  (x,y),  r'  =  (x',y' ).  Eq.  (2.17)  can  now  be  rewritten  as 
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where  is  the  field  independent  part  of  Eq.  (2.17),  i.e.,  the  right  hand  side 
with??  (t',t)  =(?7t(t',t)  =  0.  In  Eq.  (3.6)  we  exploited  the  particular  E^- 
dependence  of  Eq.  (2.17)  in  order  to  replace  -itU^E  in  the  field 

dependent  terms.  Hence  we  could  pull  the  latter  terms  out  of  the  integral  ever 
T".  No  assumption  as  to  the  smallness  of  CoV  is  involved.  Hence,  barring  a 
couple  of  differentiations,  the  problem  is  reduced  to  the  field  free  one. 

The  quantity ^ (r,r‘)  is  now  most  easily  evaluated  by  inserting  the 
integral  representations  (2.10)  and  (2.16)  and  expanding  everything  with  respect 


to  x  =  | r  -  r ' | .  The  result  is 
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One  of  the  two  sums  can  be  carried  out  in  a  multitude  of  ways,  but  the  explicit 
form  (3.7)  is  the  most  useful  one. 

We  new  have  to  insert  Eq.  (3.7)  into  Eq.  (3.6)  and  to  carry  out  the 
required  differentiations.  Needless  to  say,  Eq.  (3.6)  would  be  entirely  useless 
if  the  differential  operators  could  not  be  expanded  into  power  series  with 
respect  to We  shall  first  write  down  and  then  discuss  the  expansions  of 
exp  (-i )7l)  andvfi(ea^7);  with  the  abbreviations 


they  are 
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In  applying  the  differential  operators  (3.9)  and  (3.10)  to  Eq.  (3.7),  we  first 
note  that  each  ex'  will  effectively  be  replaced  by  a  very  small  number, 
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b  v  4  p  ^ 

-  3  r0  k°  , 


since  the  ratio  (y/E^  is  very  small  in  any  case;  for  <■%/  m  we  have  (u/E  a/  10 
for  a  laser  field  and  W/E  a,  10  for  the  radio  frequency  fields  with  J  a^  100m 


which  are  considered  in  Ref.  2.  The  effective  expansion  parameter  in  Eq.  (3.9) 
is  then 
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m  £i^ •  (3 •  10 )  with  R^  1  •  1  x  10  raixA  1  the  nuclear  radius.  Both 


quantities  are  very  small  for  experimentally  available  fields  so  that  the 
expansions  (3.9)  and  (3.10)  are  very  well  justified.  With  respect  to  Eq.  (3.7) 
we  note  that  the  main  expansion  parameter 
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is  again  small.  Hence,  it  will  be  sufficient  in  most  applications  to  keep  only 
the  lowest  power  of  x  which  is  necessary  in  view  of  the  forbiddenness  of  the 
decay.  The  sum  over  (?  is  then  restricted  to  a  few  terms. 

Let  us  first  consider  the  field  induced  enhancement  of  allowed  beta  decay. 
Since  the  expansion  parameter  (3.12)  is  in  general  smaller  than  (3.11)  we  can 
ignore  Jo(eaJlT).  Keeping  only  the  lowest  order  terms  in  Eq.(3.9)  and  using 
Bq. (3.7)  (with  x  =  0) ,  (3.6)  and  (2.8)  we  arrive  at 
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Apparently,  the  field  induced  enhancement  is  unobservably  small*.  In  Eq.  (3.14) 

we  included  the  next- to- leading  terms  in  order  to  be  able  to  make  the  comparison 

with  the  nonrelativistic  limit  of  a  fully  relativistic  calculation  [13,14]  of 

laser  enhancement  of  allowed  beta-decay.  Equation  (3.14)  agrees  with  Eq.  (113) 

of  Kef.  14.  Actually,  the  next-to-leading  terms  as  exhibited  in  Eq.  (3.14)  are 

incomplete,  because  additional  terms  resulting  from  (ea vf?)  (Eq.  (3.10))  are 

of  at  least  the  same  order  of  magnitude.  Since  the  leading  term  is  extremely 

small  already,  we  did  not  care  to  write  them  down. 

We  will  new  turn  to  a  first  order  forbidden  decay.  In  the  absence  of  the 

field  the  leading  contribution  to  the  decay  rate  oomes  from  the  term  with  k  =  1 

in  Eq.  (3.7).  Recalling  that  =  (r  -  r‘)  and  that  we  need  one  power  of  r  in 

<v/ 

each  nuclear  matrix  element  to  render  if  nonzero,  we  obtain 
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"*  previous  conclusions  to  the  contrary  [10,11]  were  flawed  and  due  to  an 
inadequate  approximation  in  evaluating  sums  over  Bessel  functions.  The 
present  authors  have  corrected  this  mistake  in  the  general  form  of  a  no-go 
theorem  [12]. 
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Considering  new  the  external  field  we  know  already  from  the  previous  example 
that  the  enhancement  with  is  formally  due  to  the  operator  TlZ,  is  extremely 
small.  Hence,  we  shall  rather  consider  the  operator  J^(ea{~z}  which,  via  its 
dependence  on  r  and  r',  generates  the  kind  of  field  induced  modifications  of  the 
selection  rules  which  are  proposed  in  Refs.  1  and  2.  Hence  in  the  presence  of 
the  field  we  do  obtain  a  non  zero  contribution  to  the  decay  rate  from  the  term 
with  k  =  0  in  Eq.  (3.7),  which  we  shall  call 
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Keeping  just  the  leading  term  in  the  expansion  (3.10)  we  get* 
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If  we  assume  that  <r_>  /<r>  =  2/3,  the  enhancement  is 
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The  last  term  in  the  square  bracket  in  Eq.  (3.10)  seems  at  first  glance  to 
yield  an  imaginary  contribution  to  the  decay  rate.  Actually  it  cancels  when 
integrated  over  with  the  nuclear  wave  functions. 


Hence  we  arrive  at  much  the  sane  conclusion  as  in  the  case  of  allowed 


beta-decay:  the  enhancement  which  is  due  to  field  induced  nod ificat ions  of  the 
selection  rules  is  by  far  too  small  to  be  observaole. 

In  order  to  obtain  the  total  enhancement  of  a  first  forbidden  decay  we 
still  have  to  add  the  modification  of  I  ,  due  to  the  operator  /jtZ ,  viz. 
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so  that  the  total  enhancement  is 


i~>  r~> 

0,0.  +  /jQ. 


) 


£?/  +  'ti  *  jr  * 

-  ^  r.)  (  r.  .  ' 


6^  D 


4-  <fta 


tC'r>  + 


(3.20) 


,  % 


To  give  a  specific  example,  let  us  look  at  the  first  order  forbidden  decay 


of  Sr  which  was  considered  in  Ref.  2.  In  that  paper  an  enhancement  of  the 


half  life  from  the  usual  29  years  to  about  10  years  is  predicted  to  occur  in  a 

3* 


radio  frequency  field  environment  with  E^8xl0  Volt/cm.  Corrected  estimates 
have  not  been  given  yet,  but  it  is  argued  [4,5]  that  enhancements  of  first  order 
forbidden  decays  should  stay  more  or  less  like  they  were.  In  contrast,  Eq. 


(3.20)  yields  for  this  case  a  relative  change  in  the  decay  rate  of  about 


2  x  10 


The  preceding  calculation  was  for  circular  polarization.  With  slightly 
more  labor,  it  can  be  done  for  linear  polarization  as  well,  with  qualitatively 
the  same  results,  i-ioreover,  it  can  be  shown  that  for  an  arbitrary  plane  wave 
field  (arbitrary  polarization  and  pulse  shape)  the  total  decay  rate  is  entirely 
independent  of  the  field  as  long  as  the  quasiclassical  limit  applies  [15]. 


IV.  CONCLUSIONS 


We  have  shown  that  the  effect  of  an  intense  plane  wave  field  cn  forbidden 
beta-decay  is  much  the  same  as  it  is  on  allowed  beta-decay.  There  are  two  types 
of  field-induced  contributions  to  the  total  decay  rate:  one  that  does,  via  its 
dependence  on  the  nuclear  coordinate,  change  the  selection  rules  of  the  decay 
and  one  that  does  not.  Both  terms  lead  to  an  enhancement  of  the  decay,  but  in 
either  case  the  enhancement  is  proportional  to  the  squared  ratio  (E/E^^  )  of 
the  applied  over  the  critical  field  E,  -  =  1.3  x  10  Volt/cm.  For  an  intense 

focused  laser  field  we  have  E/E  10~^"at  best,  for  a  radio  frequency  field 

the  ratio  is  very  much  smaller. 

Our  results  for  the  enhancement  differ  from  previous  recent  estimates 

[1. 2.4.5]  by  twenty-four  orders  of  magnitude.  The  reason  for  this  large 
discrepancy  seems  to  be  the  following:  the  original  explicit  estimates  for 
various  specific  examples  [2]  were  later  shown  to  be  faulty  [3]  and  withdrawn 

[4.5] .  However,  in  Refs.  4  and  5  the  conclusion  is  maintained  that  there  should 
be  an  enhancement  of  the  order  of  unity  without  explicit  numbers  being  given. 
This  conclusion  is  drawn  from  the  observation  that  there  is  a  significant  effect 
of  the  external  field  cn  the  differential  decay  rate.  This,  however,  does  not 
necessarily  reflect  any  substantial  changes  of  the  total  decay  rate.  We  have 
shown  previously  [12]  that  rhe  total  decay  rate  of  a  neutral  particle  decaying 
into  various  charged  particles  in  the  presence  of  an  intense  plane  wave  field, 
is  in  the  semi-classical  (WKB)  limit  entirely  independent  of  the  field.  This  is 
notwithstanding  the  fact  that  the  energy  spectrum  of  the  decay  products  is 
completely  distorted  by  the  field  as  compared  to  the  field-free  one  (  for 
examples  in  the  case  of  allowed  beta-decay  see  Ref.  13).  The  latter  effect  is 
due  to  the  essentially  classical  interaction  between  the  field  and  the  charged 
decay  products  once  the  decay  has  already  taken  place.  The  decay  itself  is  not 
affected  until  the  field  becomes  comparable  with  the  (genuinely  quantum 


E  .20 

mechanical)  critical  field  E^-^*  Reference  12  does  not  strictly  apply  to  the 
present  situation  because  the  decaying  nucleus  is  not  neutral  and  there  is  no 
allowance  for  bound  particles.  It  is,  however,  easily  extended  to  the  present 
situation  [15].  The  bottom  line  is  the  same:  enormous  effects  on  the 
differential  decay  rate  do  not  correspond  to  any  significant  enhancement  of  the 
total  rate.  The  calculation  initiated  in  Refs.  1,  2,  4,  and  5  should,  if 
finished,  lead  to  the  same  conclusion. 

We  should  add  a  critical  remark  regarding  our  calculational  approach 
exhibited  in  Section  4.  Equation  (3.6)  which  displays  the  decay  rate  as  the 
result  of  applying  a  differential  operator  to  the  free  decay  rate  is  exact  as  it 
stands.  However,  the  following  ev  luation  proceeds  via  perturbation  theory  with 
respect  to  the  external  field.  This  is  at  a  stage  where  the  quantitative 
legitimacy  of  perturbation  theory  is  obvious;  yet,  the  exact  result  is  likely  to 
contain  terms  which  are  quantitatively  minute,  but  nonanalytic  with  respect  to 
the  coupling  to  the  external  field.  This  is  known  to  be  the  case  for  allowed 
beta-decay  [13,14].  These  terms  cannot  be  obtained  by  our  procedure.  It  should 
be  emphasised  that  we  do  not  apply  perturbation  theory  to  the  differential  decay 
rate;  the  latter  containing  multiphoton  interactions  of  extremely  high  order  (up 
to  10  )  between  the  field  and  the  electron  is  clearly  inaccessible  by 
perturbation  theory  of  any  reasonable  finite  order.  However,  the  total  decay 
rate  is  a  different  matter:  for  fields  below  the  critical  field  perturbation 
theory  yields  a  safely  convergent  asymptotic  series. 

We  mention,  finally,  that  the  situation  is  entirely  different,  if  there  is 
a  significant  interaction  between  the  nucleus  and  the  field.  This  happens,  for 
instance,  if  an  intense  x-ray  laser  is  nearly  resonant  with  an  excited  nuclear 
level.  In  this  case  substantial  enhancements  of  the  total  decay  rate  are  found 
[16]. 
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A  recent  claim  that  forbidden  nuclear  beta  decay  can,  by  the  application  of  a  high-intensity  radio¬ 
frequency  field,  be  enhanced  by  many  orders  of  magnitude  is  contested.  The  effect  is  shown  to  be  nonex¬ 
istent,  at  least  within  the  theoretical  model  which  has  been  adopted  thus  far. 
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1.  INTRODUCTION 

In  two  interesting  recent  publications1'2  Reiss  claims  that 
forbidden  nuclear  beta  decay  can,  in  the  presence  of  an  in¬ 
tense,  but  readily  achievable  radio  frequency  field,  be 
enhanced  by  many  orders  of  magnitude,  the  more  so  the 
higher  forbidden  the  decay  is.  Of  course,  this  would  be  a 
fascinating  effect,  both  from  a  theoretical  and  a  practical 
point  of  view.  Unfortunately,  we  find,  at  least  in  the  frame¬ 
work  of  the  approximations  and  idealizations  underlying  this 
work,  his  conclusions  to  be  incorrect  and  the  effect  nonex¬ 
istent. 

The  approach  of  Ref.  1  to  nuclear  beta  decay  in  the  pres¬ 
ence  of  a  strong  external  electromagnetic  field  treats  the 
weak  coupling  to  first  order  of  perturbation  theory,  but  is 
supposed  to  take  into  account  the  coupling  with  the  external 
field  to  all  orders,  at  least  approximately.  Hence  in  Reiss's 
model  the  state  of  the  emitted  electron  is  described  by  the 
Volkov  wave  function,3  whereas  the  nuclear  states  are  to  be 
well  approximated  by  the  so-called  momentum  translation 
approximation  (MTA)  wave  function.4  We  shall  show  in 
Sec.  II  by  explicit  calculation  that  within  this  particular 
model  the  total  beta-decay  rate  is  essentially  independent  of 
the  external  electromagnetic  field.  Our  argument  turns  out 
to  be  independent  of  whether  or  not  we  employ  a  relativistic 
description  of  the  decay.  Hence,  for  the  sake  of  simplicity, 
we  shall  first  turn  to  a  completely  nonrelativistic  description 
within  a  long  wavelength  approximation  for  the  applied 
field.  Once  we  have  made  our  point,  the  generalization  to 
the  f_lly  relativistic  problem  is  conceptually  straightforward 
and  will  be  presented  in  Sec.  III.  Our  conclusion  that  the 
lifetime  of  a  nucleus  cannot  be  influenced  by  an  external 
field  in  the  framework  of  the  model  adopted  by  Reiss  leads 
us  to  suspect  that  Ref.  1  includes  a  calculational  error.  In 
Sec.  IV  we  point  out  such  an  error.  In  order  to  understand 


why  Reiss's  model  cannot  yield  an  enhancement  of  the  nu¬ 
clear  decay  rate  we  critically  review  the  derivation  and  justi¬ 
fication  of  the  MTA_wave  function  in  Sec.  V.  By  comparing 
the  “F*ls”  vs  “p- A”  form  of  the  interaction  Hamiltonian 
we  obtain  an  indication  of  why  the  MTA  method  is  an  un¬ 
justified  approximation.  We  then  point  out  with  the  help  of 
gauge  arguments3  that  the  MTA  wave  function  is  nothing 
but  the  correct  unperturbed  state  in  the  Coulomb  gauge.6 
This  statement  can  be  corroborated  by  writing  the  5  matrix 
for  the  beta  decay  in  different  gauges.  In  this  way  we  see 
that  in  Reiss's  theory  only  the  interaction  of  the  field  with 
the  electron  is  incorporated  exactly,  whereas  the  interaction 
with  the  nucleus  is  completely  neglected.  This  fact  is  in 
strict  contrast  to  the  statements  and  intentions  of  Ref.  1  and 
leads  to  a  completely  different  interpretation  of  the  model. 
If  the  interaction  of  the  nucleus  with  the  external  field  is 
not  incorporated  into  the  model,  a  modification  of  the  total 
decay  rate  can  only  originate  from  the  coupling  of  the  elec¬ 
tron  to  the  field.  It  was  recently  shown7  that  decays  of  neu¬ 
tral  particles  are  unaffected  by  the  application  of  optical  and, 
even  more  so,  radio  frequency  fields,  to  an  excellent  ap¬ 
proximation.  Consequently,  since  Reiss's  approach  to  the 
problem  does  not  contain  any  genuine  interaction  of  the  nu¬ 
clei  with  the  field,  it  cannot  yield  any  impact  on  the  nuclear 
lifetime.  In  Sec.  VI  we  summarize  our  various  criticisms  of 
Ref.  1. 

II.  FIELD  INDEPENDENCE  OF  THE  TOTAL 
TRANSITION  PROBABILITY 

The  starting  point  of  the  formalism  of  Refs.  1  and  2  is  the 
5  matrix  for  nuclear  beta  decay.  If  we  denote  the  weak  in¬ 
teraction,  which  causes  the  transition,  by  (gK),  the  nonrela¬ 
tivistic  limit  to  the  5-matrix  element  to  first  order  in  the 
weak  coupling  reads 


The  various  terms  in  Eq.  (2.1)  are  discussed  below.  tion  is  then  a  solution  of  the  Schrodinger  equation 

denotes  the  neutrino  wave  function,  which  is  a  plane  wave  ^  .  _ 

with  momentum  q.  Following  the  procedure  outlined  in  fe><  — lp-eA(r)]J'I'(e)(  T,t)  .  (2.2) 

Ref.  1,  we  take  the  wave  functions  of  the  charged  particles,  m 

i.e.,  the  electron  and  the  nucleus  in  the  initial  and  final  We  use  here  natural  unitsg-c-1.  Unlike  the  convention 
state,  to  be  in  the  Coulomb  gauge.  The  electron  wave  func-  in  Ref.  1  where  e  -  |e|,  we  denote  the  electric  charge  of  a 
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particle  by  e,  so  that,  for  instance,  the  electron  charge  is  e  —  —  \e  |.  The  exact  solution  of  Eq.  (2.2)  is  the  nonrela- 
tivistic  Volkov  wave  function  with  momentum  p, 


*u>(  F,/;p)-exp  -  /  £r  -  p  F  -  J  [2c A(t)p-c2Aj(t)]</t  , 


where  £  —  p1/2m,  so  that  our  nonrelativistic  theory  is  only 
consistent  for  low  energy  emitted  electrons. 

The  underlying  nuclear  model  of  Ref.  1  is  the  shell  model 
with  an  inert  0*  core,  which  is  affected  neither  by  the  beta 
decay  nor  by  the  electromagnetic  Held,  and  one  or  more 
valence  nucleons  in  an  angular  momentum  erupted  state. 
The  nuclear  wave  functions  'I'0’  and  4'</)  are  then  derived 
from  a  one-particle  Hamiltonian.  The  initial  and  final  nu¬ 
clear  states  are  approximated  in  Ref.  1  by  the  so-called 
momentum  translation  approximation  (MTA)  wave  func¬ 
tion.  This  wave  function  is  given  by 


•p’+i'm  . 


♦  F.r)-exp[/eNA(r)F]<t>0(  T.t)  ,  (2.4) 

where  eN  is  the  reduced  nuclear  charge  (Eq.  (6)  of  Ref.  11 
and  4>0  denotes  the  nuclear  wave  function  in  the  absence  of 
the  external  electromagnetic  field,  i.e., 

.  (2.5) 


V(  T)  denotes  the  nuclear  binding  potential,  and  m,  is  the 
reduced  nuclear  mass  [Eq.  (5)  of  Ref.  1).  Since  we  consider 
the  nucleus  initially  and  finally  to  be  in  an  eigenstate  <t>,(  T) 
in  —  i.f )  of  H o  with  energy  E„ ,  we  shall  later  use  d>0  in  the 
form 

<bi'"(F>r)-e"'£‘,V.(r)  .  (2.7) 

For  the  time  being  we  shall  adopt  the  MTA  wave  function 
in  the  same  way,  as  it  is  used  in  Ref.  1,  and  postpone  a  de¬ 
tailed  discussion  on  this  subject  to  Sec.  V. 

The  total  transition  probability  per  unit  time  is  calculated 
from  the  5 -matrix  element  (2.1)  by 


If  we  insert  the  5-matrix  element  (2.1)  into  Eq.  (2.8)  using 
the  wave  functions  (2.3),  (2.4),  (2.7),  and  plane  waves  for 
the  neutrino,  the  total  rate  f  takes  the  form 


f  -  ^  t  jljm^  y  J*  d3r  d V  J_  fn  dt  df'expl  /£0(  /'  - t )  ]  exp  ( -  ie  ( A  ( t‘ )  F'  -  A  ( / )  T 1 1 


x  G,„,(  F'.r';  F,r  )G<e>{  F',r';  T.t ) 


x  \4>}i  DigVHA  F)H*,*(  F')(glO*<M  F’)l  . 


Here  £o —  £<-£/ denotes  the  nuclear  energy  released  in  the  beta  decay.  We  furthermore  used  the  relation  e(-e/-e  [see 
Ref.  1,  Eq.  (49),  and  our  convention  for  the  sign  of  e].  In  fact,  this  relation  is  only  approximate,  since  the  minor  impact  of 
the  external  Held  on  the  nuclear  core  is  neglected.  The  sum  over  the  neutrino  states  is  expressed  in  terms  of  a  plane  wave 
Green’s  function  with  a  dispersion  £,-  |q|, 

G(,»(  F',»';F,r)-  fd3q  ¥*,,(  F,r;q )¥(,)(  F',r';q)-  J dyq  exp[  -  »  |q|(r'- 1 ))exp(/q(  F'-  F)J  .  (2.10) 

All  the  electronic  contributions  to  the  total  rate  V  are  contained  in  the  nonrelativistic  Volkov  Green’s  function 
G(«»(F\r';F.r)-  § d'pVM{  F.r;p)'*r(e»(  F’,r';p) 

"eXP|''^^  XJ(r)rfr)/^CXP(_'[i2^ipl-(r'“r+m  X  X(rHP|  '  (2  ,,) 

The  integral  (2.11)  can  be  reduced  to  a  Gaussian  integral  and  we  find 


[  2m v  I  __(m(F'-F):|  F'-F  f «' t /  v,l 
- -  exp  /- - ; -  exp  ie — - -  I  A(r)dr 

l»(r  -/)  J  [2  t  -t\  t  -t  J<  J 

*  exp(  - ^  -  ^7~(  A(r)rfr)  |  . 


The  AJ  term  in  Eq.  (2.11)  contributes  to  an  effective  mass  (see  Sec.  III). 

The  essential  fide  dependence  of  the  total  rate  T  is  concentrated  in  the  factors 

exp(*[A<i')F’-A(l)Fl| 
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and  the  Volkov  Green's  function  G <,).  In  Ref.  1  the  entire 
effect  of  changing  the  degree  of  forbiddenness  of  the  nu¬ 
clear  decay  is  derived  from  the  factors  exp< ieA7)  in  Eq. 
(2.4).  We  shall  now  show  explicitly  that  these  factors  actu¬ 
ally  cancel  out  of  the  total  decay  rate  against  corresponding 
terms  in  G(e>  and  hence  cannot  give  rise  to  the  effects  cal¬ 
culated  in  Ref.  1 . 

By  inserting  Eqs.  (2.12)  and  (2.10)  into  the  total  decay 
rate  (2.9),  we  see  that  in  the  field-free  limit  the  major  con¬ 
tribution  to  the  time  integral  in  Eq.  (2.9)  comes  from  time 
differences  t'—t ,  for  which  the  phase  of  the  integrand  is  sta¬ 
tionary: 


If'-rl- 


w(r-r)2  1/2 
2(£0-£.) 


(2.13) 


The  estimate  (2.13)  is  not  significantly  affected  by  the 
external  field,  in  particular,  when  there  is  no  energy  transfer 
from  the  field  to  the  electron,  as  Reiss  assumes. 

It  was  shown  in  Ref.  7  that  in  the  classical  limit  K  —  0 
only  times  r'  infinitesimally  close  to  t  contribute  to  the  total 
transition  rate  T.  In  this  limit  the  field-dependent  exponen¬ 
tials  in  the  integrand  of  Eq.  (2.9)  are  unity  and  only  the 
field-independent  part  survives.  Thus  in  the  classical  limit 
the  total  decay  rate  would  be  unaffected  by  the  external 
field.  This  is  in  agreement  with  the  results  of  Ref.  7. 

The  remaining  field  dependence  in  T  is  then  only  due  to 
quantum  effects.  It  will  be  found  to  be  very  small  for  the 
parameters  of  Ref.  1,  as  we  shall  now  demonstrate.  We  can 
rewrite  the  space-  and  field-dependent  exponential  in  the 
electron  Green’s  function  (2.12)  with  the  help  of  integration 
by  parts  and  find 


exp 


ie—^ - — ~  dr  A(t)  -  exp|/e(A(r')  F'  — A(r )  r]!exp[/e  dT  R(t  )E(r)j  , 


with 


(2.14) 


R(r)-F'--i-7 — — (f'-r)  .  (2.15) 

t  -t 

By  inserting  Eq.  (2.12)  with  Eq.  (2.14)  into  Eq.  (2.9),  we  realize  that  the  phase  factors  exp(;eA  F)  cancel  in  the  total  decay 
rate  and  we  obtain 

r  -  ::  lim  \  (d}r  d}r'  Cdt  dt'expliE0(t'-l)]G{,)(T'.l'-,T.l) 

(2jt)®  r~m  T  J  J 

exp(/e  J'  dr  R(t )E(r ) J 

XeXP(_'^[X'  drAHr)--^ [£  drMr)]  ]| 

x[d>;(F)(gl')</.,(F)](«;(F')(gK)*<f./(F’)]  .  (2.16) 


exp/f  If piy 

i(l-t)  2  t'-t 


In  Ref.  1,  Eq.  (74),  the  nuclear  intensity  parameter  r. 


z  -  (elAlRo)1  , 


(2.17) 


which  specifies  the  magnitude  of  the  phase  expOe  A  F),  was 
assumed  to  be  of  order  of  unity.  The  actual  remaining  ex¬ 
ponential  with  a  field  and  space  dependence  in  Eq.  (2.16) 
contains  the  integral 


1 1  J^'  R(r)E(r)</r|~eR0lElU'- 

—  eR  o«u  I A 


mR  o 


E  o 


1/2 


(2.18) 


Here  R0  denotes  the  nuclear  radius,  and  the  length  of  the 
time  interval  was  estimated  by  Eq.  (2.13).  We  furthermore 
determined  the  amplitude  of  the  electric  field  by  I E I  —  &»  |  A I , 
corresponding  to  a  monochromatic  plane  wave  with  fre¬ 
quency  w.  and  applied  Eq.  (82)  of  Ref.  1.  Equation  (2.18) 
shows  that  the  field-  and  space-dependent  exponential  in 
Eq.  (2.16)  is  unity  to  an  excellent  approximation,  i.e.,  the 
factors  exp(ieAF),  which  are  the  origin  of  the  large 
enhancement  obtained  in  Refs.  1  and  2,  cancel.  The  A 1 


I - — - 

terms  in  Eq.  (2.16)  do  not  depend  on  F  and  can  therefore 
not  change  the  degree  of  forbiddenness  of  the  decay. 

We  also  see  from  the  estimate  (2.18)  that  the  actual 
parameter  that  governs  the  field  impact  on  the  nuclear  life¬ 
time  is 

z'- (e<u|A|/?o  ):- z(uiRo)1  ,  (2.19) 

which  is  much  smaller  than  z.  This  means  that  one  needs  a 
much  higher  field  intensity,  or  fields  with  a  much  shorter 
wavelength  than  the  one  applied  in  Refs.  1  and  2,  to  pro¬ 
duce  a  noticeable  effect  of  the  external  field  on  the  nuclear 
lifetime. 


III.  RELATIVISTIC  THEORY 

We  shall  now  address  ourselves  to  a  relativistic  treatment 
of  the  electron  in  analogy  with  Ref.  1.  For  a  quantitative 
analysis  this  is  indispensable  because  the  electron  energy 
will,  in  general,  be  relativistic,  and  because  the  entire  decay 
process  is  intrinsically  relativistic.  However,  the  crucial 
point  of  the  preceding  analysis,  the  actual  replacement  of 
the  superficial  appearing  gauge  factors  exp(ieAF)  by 
explie  fdr  R(r)E(r)]  in  the  total  decay  rate,  will  proceed 
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completely  analogously  in  the  relativistic  treatment.  We 
then  shall  have  to  examine  the  remaining  expressions  for 
possible  additional  terms  incorporating  AF  which  might 
yield  the  large  enhancements  claimed  by  Reiss.  These  will 
have  to  be  relativistic  quantum  terms  induced  by  a  radio 
frequency  field.  Hence  it  is  not  likely  that  they  will  play  any 
significant  role,  and,  in  fact,  they  will  not  turn  out  to  do  so. 

The  starting  point  will  now  be  the  relativistic  S -matrix 
element  Sj i  as  given  in  Eq.  (7)  of  Ref.  1.  As  we  did  in  Eq. 
(2.8),  we  shall  again  sum  over  the  electron’s  final  momenta 
p  and  spins  s,  obtaining 


-  f  d*pO{p0)S(p*-m1)  JlSyil2 

fd*xd*x' l^/(x)yM(l-Ky5)'l'/U)] 

4m  ** 

x(¥,(x')(l  +xy,)y^/(x')] 
x[*<,,(x')(l  +  y  s)y°G(*\x‘,x) 

xr»‘(l-y5)^(,)(x)l  .  (3.1) 

Here 

G{  *  Hx'.x )  -  J d*p  0(po)8(p2—  m2)E (x'.p )  (jf  +  m  )E ( x.p ) 

(3.2) 

is  the  homogeneous  positive  frequency  Volkov  Green’s 
function  replacing  Eq.  (2.11),  and  we  have  written  the  Vol¬ 
kov  solution  in  the  form 

'V(,Ax,p)-  E(x.p)u(p.s)  , 

where  u  (p,s )  is  the  free  Dirac  spinor  so  that 

(p- m)u(p.s)  —  0  . 

We  are  again  using  the  Coulomb  gauge  so  that  the  nuclear 
wave  functions  are  given  by  Eq.  (2.4). 

The  quantity  R,  when  summed  over  the  final  states  of  the 
nucleus  and  the  neutrino,  yields  the  total  decay  rate  of  Eq. 


(2.8).  The  decisive  point  is  that  any  Green’s  function  in  the 
presence  of  an  external  field  AM(x)  can  be  split  as* 

Gn(x\x)-4>(x',x)Gn(x'.x)  ,  (3.3) 

where 

0(x',x)-exp[-/e  dx^AUx)}  (3.4) 

is  Schwinger’s  line  integral  which  is  to  be  integrated  over  a 
straight  line  connecting  x  with  x  (the  definition  can  be 
given  in  a  path-independent  way  (see  Ref.  8)],  and 
G(  ’(x'.x )  is  gauge  invariant,  i.e.,  depends  only  on  E  and  S. 
Hence  the  entire  gauge  dependence  is  isolated  in  the  line  in¬ 
tegral  (3.4).  Collecting  now  all  the  fA-dependent  ex¬ 
ponentials  in  Eq.  (3.1)  we  obtain 

exp(/e[  FA(/)- F'A(/')l-te  dXpAHx)\  . 

(3.5) 

Evaluating  the  line  integral  yields 

f*  dxMAM(x)  —  —  — -p7 — ffa  d(kx)A(kx)  , 

J*  k(x  —x)  Jkx 

(3.6) 

where  denotes  the  wave  vector  of  the  field  AM(kx).  In 
the  long  wavelength  approximation,  which  is  finally  also 
adopted  in  the  relativistic  treatment  of  Ref.  1,  Eq.  (3.6) 
reproduces  the  expression  which  we  already  encountered  in 
the  nonrelativistic  treatment.  The  exponential  (3.S)  is  then 
identical  to  the  corresponding  exponential  in  Eq.  (2.9)  with 
the  nonrelativistic  Green’s  function  (2.12),  and  the  argu¬ 
ments  pointed  out,  following  Eq.  (2.12),  apply  to  the  rela¬ 
tivistic  case  as  well. 

Finally,  in  order  to  make  sure  that  the  gauge-invariant 
remainder  G<  +  >  in  Eq.  (3.3)  does  not,  so  to  speak  by  the 
back  door,  reintroduce  corresponding  exponentials,  we  shall 
now  write  down  the  complete  Green’s  function  G'^M x'.x). 
The  easiest  way  to  obtain  it  is  by  straightforward  evaluation 
of  Eq.  (3.2),  given  the  Volkov  wave  functions  E(x,p).  The 
analogous  approach  has  been  followed,  e.g.,  in  Ref.  9,  in 
the  case  of  the  Feynman  propagator.  The  result  (for  arbi¬ 
trary  polarization,  i.e..  A*- 2}mlefA,l( ).  (mmkx,  Ace,- 0, 
e,ey-  - 6(/)  is  very  closely  related  to  the  former  and  reads 


<j‘*’(x\*)-  f"  —*(m2+n-/—  ,*-)2 

J --  sl  4j 


x{9[s(f'-f)l  + 


-  /,((' -OLt  +  i  <{' -  f  )y  +  H 


-(x"-x')L,  +  s(f-{)U,I-A/,J)--^ 


sT 


+  »y5<'[-«(/U"-x')A/y  +  2j(f'-f)ttfi,A/y]J  +  /sgn(s)/8({'-f)|  .  (3.7) 


Here  x‘  denotes  the  two  vector  components  of  x“  transverse  to  k We  adopt  a  summation  convention  for  the  indices  »  and 
J  extending  from  1  to  2,  and  -«*,  ««-  1.  The  functions  Ls  A*,  u.td  T  are  given  by 


1128 


COMMENTS 


29 


Jf  J*  d-n'A,(n') 

2«w)  hu,u(<)-?^?  -f dvMv)\  ■ 


(3.8) 

(3.9) 
(3.10) 


F.5  _ 


We  incidentally  note  that  the  quantity  T  ((',()  in  view  of 
Eq.  (3.7)  specifies  a  space-time-dependent  mass  correction. 
It  already  showed  up  in  the  nonrelativistic  treatment  [cf.  the 
last  exponential  in  Eq.  (2.12)]. 

We  now  have  to  check  whether  Eq.  (3.7)  contains  any 
terms  similar  to  cos(eTF-ea  F'),  which  is  the  sole  cause 
of  the  enhancement  in  Ref.  1  (cf.  Eq.  (110)].  There  are 
two  candidates  in  Eq.  (3.7): 

-X(x"-x')Li.  -iysKtyW-x^Mj  . 

We  notice  that  they  are  both  only  linear  in  (x"-x‘). 
Hence  they  can  at  best  reduce  the  order  of  forbiddenness  by 
one.  Since  they  are  proportional  to  X,  their  order  of  magni¬ 
tude  is  <ue|A|7?o“  which  has  to  be  compared  with  the 
electron  mass  m  in  the  square  bracket  in  Eq.  (3.7).  Since 
a>/m  «  1,  in  particular,  fora  radio  frequency  field,  whatev¬ 
er  effect  is  caused  by  these  terms  will  be  very  small,  as  was 
to  be  expected. 

IV.  AN  ALGEBRAICAL  ERROR  IN  REF.  1 

The  above  considerations  strongly  suggest  that  there  is  a 
computational  error  in  Ref.  1.  Hence  we  have  cursorily 
checked  some  of  the  calculations.  In  fact,  it  appears  that 
factors  ( -i)J  and  /"  are  missing  on  the  right-hand  side  of 
Eqs.  (56)  and  (57)  of  Ref.  1,  respectively,  which  should 
read 

exp(  -  ieA  F)  -exp[  -  tea  F cos (kx  +p)] 

-  J  expl  —  ijikx  +p)](  —  i)J 

(4.1a) 

and 

exp(/eA'F')  — exp(/ea  F'cos(Aar'+p)J 

—  5)  y.(ea F')exp[//?t (Aar’+p)]/"  .  (4.1b) 

These  factors  seem  to  be  consistently  missing.  If  they  are 
included,  the  crucial  equation  (110)  of  Ref.  1, 


I - 

V.  REVIEW  AND  CRITICISM  OF  THE  MTA 
WAVE  FUNCTION 

For  a  review  of  the  MTA  wave  function  from  Reiss’s 
point  of  view  we  refer  to  Ref.  4.  The  MTA  method  claims 
to  be  “gauge  specific’’  for  the  Coulomb  gauge  (C  gauge; 
vector  potential  A  and  vanishing  scalar  potential  A o— 0). 
The  Schrodinger  equation  in  the  long-wavelength  approxi¬ 
mation  for  the  vector  potential  A  then  has  the  form 

i-f ¥(F.r)-[//o+//10)WF,r)  ,  (5.1) 

ot 

with 

//o”  v~P1+  V(  F),  H\(t) - -A(r)p+-^-AJ(r)  . 

2m  m  2m 

(5.2) 

If  one  writes  an  ansatz  for  ¥(  F,r)  in  the  form 

¥(  F,r)  «exp[/eA(f)  F]<t>(  F.t)  ,  (5.3) 

then  the  new  wave  function  <M  F.r )  obeys  the  equation  of 
motion 

/-f-<J>(  F.t)- [//o  +  /fj(i)N>(F.r)  .  (5.4) 

0/ 

with 

//,(!)- -eE(r)F  .  (5.5) 

The  so-called  “momentum  translation  approximation” 
consists  of  neglecting  the  perturbation  Hi  in  Eq.  (5.4)  and 
replacing  <!>  in  Eq.  (5.3)  by  the  unperturbed  wave  function 
4>0,  given  by 

F.t )-//0<!>o(F,r)  .  (5.6) 

Ot 

We  then  obtain  the  MTA  wave  function  as  an  approxima¬ 
tion  of  the  exact  solution 

¥<>(  F.t ) -exp[ieA(f )  F]4>0(  F.f )  .  (5.7) 


2<-)V„(eaF-eaF')-cos(eaF-faF')  ,  (4.2) 

« 

loses  the  factor  (  -  )f  on  the  left-hand  side  and  is  changed 
into 

Jy^fea  F-ea  F')- 1  .  (4.3) 

« 

and  there  is  no  field  induced  enhancement  of  forbidden 
beta  decay,  in  agreement  with  the  arguments  previously  put 
forward  in  this  paper. 


This  procedure  is  considered  to  be  justified  if  H2  is  much 
smaller  than  Hi,  in  particular,  when  Hi  is  too  large  to  be 
treated  as  a  perturbation  with  respect  to  Ho,  whereas  H j  still 
is  a  small  perturbation  compared  with  Ho  (i-e.,  the  magni¬ 
tude  of  H 2  is  small  as  compared  with  a  characteristic  energy 
of  the  field-free  problem).10  The  condition  "Hi  small  com¬ 
pared  with  Hi"  means  that  the  transition  matrix  elements 
from  an  initial  unperturbed  state  |i)  to  a  final  unperturbed 
state  I/)  are  much  smaller  when  the  transition  is  induced 
by  the  residual  interaction  Hi  instead  of  Hi.  The  states  It) 
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and  |/>  are  eigenstates  of  Ho- 


Ho\n  >-£„|n>  . 

(5.8) 

By  using  the  commutator  relation 

m 

(5.9) 

as  well  as 

lE|-w|A|  . 

(5.10) 

we  find 

K/ltfJ(r)|,)|  0,  • 

|</|//,(r)|,)|  E,-Ej  ‘ 

(5.11) 

Hence  in  situations  where  the  field  frequency  is  much 
smaller  than  the  considered  transition  energy,  the  MTA 
procedure  might  seem  to  be  well  justified. 

As  we  have  noted,  central  to  the  MTA  is  the  condition 
that  the  transition  matrix  elements  for  the  interaction  Hi 
are  much  smaller  than  for  its  counterpart  Ht.  This  applies, 
in  particular,  to  nuclear  transitions  in  the  presence  of  an  op¬ 
tical  or  even  radio  frequency  field.  But  one  must  be  careful 
in  drawing  conclusions  from  or  making  approximations 
based  upon  the  estimate  <5.1 1).  Instead  of  transition  matrix 
elements,  one  actually  has  to  compare  transition  probabili¬ 
ties,  which  are  the  directly  measurable  quantities.  As  we 
shall  see  below,  the  factor  a >HE,-Ej)  does  not  occur  after 
integration  over  time  and  the  two  interactions  H\  and  Hi 
give  the  same  transition  probabilities  for  any  ratio  of  <u  over 
(£,-£>). 

Let  us  repeat  a  supposedly  well  known  argument.  In  the 
interaction  picture  the  transition  amplitude  to  first  order  of 
perturbation  theory  reads 

F(  T)  -  -  /  f'di  expl  -  HE,  -  Ef)t  ]  (f\H,M(t)\i)  . 

(5.12) 

If  we  use  the  interaction  H,ni“Hi,  we  find 

F(T)-/e(/|r|/>  J]  dr  expl  - /(£,-£,)/ )!(/)  ,  (5.13) 

For  the  case  H,nl-H i,  let  us  consider  the  situation  where 
the  field  is  switched  on  at  time  t  -  0  and  that  we  are  only 
interested  in  the  transition  probability  at  time  T,  when  the 
field  is  already  switched  off,  so  that 

A(0)-0.  A(T)-0  .  (5.14) 

This  is  the  usual  experimental  situation.  We  can  then 
rewrite  the  integral  in  Eq.  (5.12)  by  partial  integration 

Ja  exp[  -  »'(£,- E/)t]A(t)dt 

~  J*  exp[  -/(£,  -E/)t  ]E(t)dt  .  (5.15) 

L,  ~  L/ 

We  have  now  obtained  a  factor  (£,  -  E/)  by  replacing  A  by 
I  in  the  integrand,  instead  of  the  factor  u,  which  entered 
the  relation  (5.11)  via  Eq.  (5.10).  By  using  Eqs.  (5.9)  and 
(5. IS)  we  find  for  //,„,-// 1  again  the  transition  amplitude 
(5.13). 

This  simple  calculation  shows  that  the  transition  probabili¬ 
ty  is  the  same  whether  one  uses  the  interaction  H\  or  Hi. 


This  holds  true  to  any  order  in  the  perturbation  series.11 
We  can  furthermore  drop  the  restriction  (5.14)  and  obtain 
equal  transition  amplitudes  F(T)  at  any  lime  T,  if  we  take 
into  account  that  the  same  physical  state  is  represented  by 
different  state  vectors  in  different  gauges  if  5(0)* 0  or 
A(T)*0  (see  next  paragraphs).  The  error  will  therefore 
be  of  the  same  order  of  magnitude  whether  one  neglects  the 
interaction  Hi  in  Eq.  (5.4)  or  H\  in  Eq.  (5.1)  for  any  ratio 
of  the  field  frequency  <a  over  the  typical  transition  energy 
E,-E/.  Since  the  approximate  solution  <J>0  of  Eq.  (5.4)  is 
an  unperturbed  wave  function,  these  arguments  indicate 
that  also  the  MTA  wave  function  represents  a  noninteract¬ 
ing  slate. 

In  order  to  see  that  the  MTA  wave  function  represents 
nothing  but  a  noninteracting  state,  we  must  consider  the 
different  gauges  involved  in  the  problem.  The  wave  func¬ 
tions  of  all  charged  particles  participating  in  any  reaction  or 
decay  process  must  be  taken  consistently  in  the  same  gauge. 
The  Coulomb  gauge  is  convenient  for  the  calculation  of  the 
Volkov  wave  function  of  a  free  particle  in  the  presence  of 
an  external  electromagnetic  field.  To  derive  the  properties 
of  the  nonrelativistic  MTA  wave  function  it  is  more  instruc¬ 
tive  to  begin  with  the.  electric  field  gauge  [E  gauge;  vanish¬ 
ing  vector  potential  A-0  and  scalar  potential  -eE(/)FJ. 
Wave  functions  are  transformed  from  the  C  gauge  (notation 
♦ )  to  the  E  gauge  (notation  <P )  by  the  unitary  transforma¬ 
tion 

<b(  F,r)-  expl  -ieA(r)  r'l'H  f*,r)  ,  (5.16) 

where  A  is  the  vector  potential  in  long  wavelength  approxi¬ 
mation  for  the  C  gauge. 

The  Schrodinger  equation  in  the  E  gauge  has  the  form 
(5.4),  which  was  derived  at  the  beginning  of  this  chapter  in 
a  different  context.  If  we  entirely  neglect  the  interaction  of 
the  particle  with  the  field  to  a  zeroth  approximation,  we  ob¬ 
tain  from  Eq.  (5.6)  as  an  approximate  solution  the  wave 
function  <J>0  in  the  absence  of  the  external  field. 

In  the  present  example  of  beta  decay,  the  nuclear  state 
has  to  be  determined  from  the  Schrodinger  equation  in  the 
C  gauge,  i.e.,  from  Eq.  (5.1),  since  the  electron  wave  func¬ 
tion  is  given  in  the  C  gauge.  If  we  approximate  the  solution 
of  Eq.  (5.1)  by  the  noninteracting  state  <l>o  in  the  E  gauge 
and  use  the  gauge  transformation  (5.16),  we  obtain  the 
wave  function  'l'0  of  Eq.  (5.7),  the  MTA  wave  Junction. 
Although  ¥0  depends  on  the  vector  potential  A,  it  still 
represents  a  noninteracting  state  since  a  state  vector  which 
does  not  include  any  interaction  with  the  external  field  in  a 
particular  gauge  (here  the  E  gauge),  neither  does  so  when 
transformed  to  any  other  gauge  (here  the  C  gauge). 

It  should  be  mentioned  that  the  wave  function  <J>0  is  not 
the  correct  noninteracting  solution  in  the  C  gauge.  In  the  C 
gauge  the  operator  of  the  canonical  momentum  p  -  —  i  ^ 
differs  from  the  operator  of  the  kinetical  momentum. 
Hence  in  the  C  gauge  the  Hamiltonian  Ho  in  Eq.  (5.2)  does 
not  describe  the  situation,  where  the  interaction  of  the  parti¬ 
cle  with  the  field  is  entirely  neglected.  The  correct  pro¬ 
cedure  is  first  to  transform  from  the  C  gauge  via  (5.16)  to 
the  E  gauge,  where  the  vector  potential  vanishes  and  the 
operators  of  the  kinetical  and  canonical  momentum  are 
identical,  so  that  p1/2m  is  the  operator  of  the  kinetical  ener¬ 
gy.  Hence  in  this  gauge  the  Hamiltonian  in  Eq.  (5.6)  really 
specifies  the  field-free  motion,  and  4>o  represents  the  state 
in  the  absence  of  the  field.  The  wave  function  in  the  C 
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gauge,  which  completely  neglects  the  interaction  with  the 
field,  is  then  given  by  Eq.  (5.7),  i.e.,  by  the  MTA  wave 
function.  For  a  lucid  and  thorough  discussion  of  gauge  in¬ 
variance  in  quantum  mechanics  see  Ref.  5. 

Applying  this  discussion  to  the  nuclear  beta  decay,  we  see 
that  the  S  matrix  (2.1)  includes  the  exact  electron  state  and 
noninteracting  nuclear  states  in  the  C  gauge.  In  the  nonrel- 
ativistic  approach  the  S  matrix  could  also  be  written  in  the  E 
gauge.  The  unperturbed  nuclear  state  is  then  given  by  the 
field-free  wave  function  <J>o,  and  the  gauge  phase 
expGeAF)  is  incorporated  in  the  corresponding  Volkov 
wave  function,  which  we  obtain  from  Eq.  (2.3)  and  the 
transformation  (5.16).  The  entire  field  dependence  is  then 
concentrated  in  the  electronic  wave  function.  This  pro¬ 
cedure  yields  the  same  analytical  form  of  the  S -matrix  ele¬ 
ment  as  Eq.  (2.1). 

The  factor  exp(ieAF)  always  appears  when  one  con¬ 
sistently  combines  the  nuclear  state  in  the  absence  of  the 
field,  which  is  related  to  the  E  gauge,  with  a  Volkov  state, 
which  can  be  calculated  in  a  simple  way  for  the  C  gauge. 
This  factor  can  be  shifted  from  the  nuclear  to  the  electronic 
wave  function  in  different  gauges  and  has  nothing  to  do 
with  the  interaction  of  the  nucleus  with  the  field.  The  iden¬ 
tification  of  the  factor  exp(ieAF)  in  Eq.  (5.3)  as  a  gauge 
transformation  was  in  the  present  context  first  pointed  out 
in  Ref.  6.  In  Ref.  4,  Reiss  does  not  consider  this  factor  as  a 
gauge  transformation  but  rather  calls  it  a  “unitary  transfor¬ 
mation  within  the  Coulomb  gauge.”  He  then  claims  that 
due  to  this  phase  factor  the  MTA  wave  function  fairly 
represents  the  effects  of  the  applied  field  on  the  particle  to 
any  order  of  interaction.  Relying  on  this  interpretation, 
Reiss  attempts  in  Ref.  1  to  derive  the  entire  effect  of  chang¬ 
ing  the  degree  of  forbiddenness  of  a  nuclear  decay  and  of 
enhancing  nuclear  decay  rates  from  this  factor  exp(ieAF). 

If  the  interaction  of  the  nucleus  with  the  field  is  neglect¬ 
ed,  only  the  coupling  of  the  electron  with  the  field  can  still 
modify  the  nuclear  lifetime.  But  as  it  was  shown  in  Ref.  7, 
this  can  only  happen  for  a  field  which  is  much  stronger  than 
the  one  considered  in  Ref.  1,  as  long  as  the  field  frequency 
is  very  small  as  compared  with  the  nuclear  decay  energy. 
Therefore  the  model  used  in  Ref.  1  cannot  result  in  an  ap¬ 
preciable  change  of  the  nuclear  lifetime. 


V.  SUMMARY 

The  basic  formalism  developed  at  length  in  Ref.  1  is 
essentially  the  same  as  in  Refs.  12  and  13  (cf.  Ref.  14). 
The  theory  of  Ref.  1  only  differs  from  this  previous  work  by 
using  the  correct  noninteracting  wave  function  'Po  (2.4)  in¬ 
stead  of  d>o  (2.5)  (but  with  the  misguided  intention  of  incor¬ 
porating  the  interaction  with  the  field)  and  by  considering 
linear  polarization  of  the  field  instead  of  circular  polariza¬ 
tion.  Whether  the  field  is  linearly  or  circularly  polarized  is 
rather  immaterial  for  the  effects  in  question.  Considering 
linear  in  place  of  circular  polarization  mainly  increases  the 
calculational  labor  without  adding  additional  insight. 

We  would  like  to  concentrate  on  two  central  objections  to 
Refs.  1  and  2:  (i)  Contrary  to  its  intentions,  Reiss’s  model 


does  not  include  any  interaction  between  the  nucleus  and 
the  field.  Notably,  the  so-called  momentum-translation 
(MTA)  wave  function  is  nothing  but  the  correct  free  wave 
function  in  the  Coulomb  gauge,  as  was  shown  in  Sec.  V  in 
three  different  ways,  (ii)  The  total  decay  rate  cannot  be  in¬ 
fluenced  by  an  external  field  as  chosen  in  Ref.  1,  if  one 
describes  the  electron  by  the  relativistic  Volkov  wave  func¬ 
tion  and  the  nucleus  by  the  nonrelativistic  MTA  wave  func¬ 
tion.  This  can  be  derived  by  general  arguments  from  point 
(i)  and  Ref.  7.  We  also  prove  this  explicitly  in  Secs.  II  and 
III.  There  we_show  that  in  a  correct  treatment  the  gauge 
factors  exp(/eA  F)  (from  which  Ref.  1  derives  its  entire  ef¬ 
fect)  mutually  cancel  in  the  total  decay  rate.  Our  argument 
shows  that  this  cancellation  is  independent  of  the  polariza¬ 
tion  and  pulse  shape  of  the  external  field.  The  dramatic 
enhancements  of  Refs.  1  and  2  seem  to  be  due  to  an  alge¬ 
braic  error.  When  this  is  corrected  no  enhancement 
remains,  in  agreement  with  the  results  of  Secs.  II  and  III. 

The  physical  concept  underlying  Refs.  1  and  2  differs 
from  the  previous  work  (Refs.  12  and  13)  by  intending  to 
concentrate  on  forbidden  beta  decay  in  the  presence  of  an 
intense  radio  frequency  field  rather  than  an  optical  field. 
Some  final  remarks  on  these  two  new  aspects:  (a)  Chang¬ 
ing  the  degree  of  forbiddenness  of  a  nuclear  decay  by  modi¬ 
fying  the  nuclear  states  under  the  impact  of  an  external  field 
is  an  exciting  idea,  but  cannot  be  achieved  by  Reiss's 
model,  which  only  uses  noninteracting  nuclear  states.  In 
particular,  it  should  be  emphasized  that  Reiss  tries  to  treat 
forbidden  beta  decay  along  the  same  lines  as  allowed  beta  de¬ 
cay,  i.e.,  a  multipole  expansion  of  the  lepton  wave  functions 
and  relativistic  corrections  to  the  nuclear  wave  functions, 
both  of  which  normally  enable  forbidden  decays  to  take 
place,  are  not  considered,  since  the  dominant  mechanism  of 
the  enhancement  is  supposed  to  originate  from  the  gauge 
factors  exp(ieAF).  (b)  The  advantage  of  radio  frequencies 
as  compared  with  optical  frequencies  seems  to  lie  in  the  fact 
that  larger  values  of  the  parameter  (ea/m  )2  can  be  achieved 
at  much  lower  field  strengths.  But  if  very  high  intensity  ra¬ 
dio  frequency  fields  with  a  wavelength  of  X  — 100  m  (Ref. 
2)  are  applied,  the  applicability  of  the  Volkov  solution, 
which  assumes  a  plane  wave  field  of  infinite  extent  in  space 
and  time,  seems  very  doubtful  and  certainly  requires  some 
justification.  A  relativistic  electron  which  moves  freely  in 
such  a  field  performs  an  oscillatory  motion  over  a  distance 
of  X.  Such  electromagnetic  fields  also  raise  experimental 
problems,  since  the  atomic  electrons  tend  to  shield  an  exter¬ 
nal  low  frequency  field,  reducing  its  field  strength  at  the  site 
of  the  nucleus  by  orders  of  magnitude.15 

We  would  finally  tike  to  emphasize  that  it  remains  an 
open  question  as  to  whether  properly  including  the  interac¬ 
tion  between  the  nucleus  and  the  field  might  yet  lead  to 
some  enhancement  of  forbidden  beta  decay,  although  we 
believe  that  the  orders  of  magnitude  stated  in  Refs.  1  and  2 
are  very  unlikely  to  be  achieved  by  the  latter  mechanism. 

We  benefited  from  discussions  with  D.  H.  Kobe  and  P. 
Zoller.  This  work  was  supported  in  part  by  the  Air  Force 
Office  of  Scientific  Research  under  Contract  No.  AFOSR- 
81-0128. 
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